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SUMMARY 


Previous work at NASA Langley Research Center (LaRC) involved fabrication and 
testing of composite beams with embedded, pre-strained shape memory alloy (SMA) ribbons 
within the beam structures. That study also provided comparison of experimental results with 
numerical predictions from a research code making use of a new thennoelastic model for shape 
memory alloy hybrid composite (SMAHC) structures. The previous work showed qualitative 
validation of the numerical model. However, deficiencies in the experimental-numerical 
correlation were noted and hypotheses for the discrepancies were given for further 
investigation. The goal of this work is to refine the experimental measurement and numerical 
modeling approaches in order to better understand the discrepancies, improve the correlation 
between prediction and measurement, and provide rigorous quantitative validation of the 
numerical analy sis/ design tool. The experimental investigation is refined by a more thorough 
test procedure and incorporation of higher fidelity measurements such as infrared 
thermography and projection moire interferometry. The numerical results are produced by a 
recently commercialized version of the constitutive model as implemented in ABAQUS and 
are refined by incorporation of additional measured parameters such as geometric imperfection. 
Thennal buckling, post-buckling, and random responses to thermal and inertial (base 
acceleration) loads are studied. The results demonstrate the effectiveness of SMAHC 
structures in controlling static and dynamic responses by adaptive stiffening. Excellent 
agreement is achieved between the predicted and measured results of the static and dynamic 
thermomechanical response, thereby providing quantitative validation of the numerical tool. 



Chapter 1 


INTRODUCTION 


1.1 Shape Memory Alloys 

Shape Memory Alloys (SMA) have become increasingly researched over the past 30 
years, due to their fascinating functional properties. In the 1950’s the Shape Memory Effect 
(SME), the recovery of defonnation to an original shape by heating to a critical temperature, 
was discovered in Gold (Au 47.5%) and Cadmium (Cd 52.5%) Alloy by Chang and Read 1 . 
Then in 1963 the SME was discovered at the Naval Ordinance Laboratory in equiatomic 
composition Titanium-Nickel Alloy (Ni 52%-57%) by Beuhler and Wiley (Nitinol) 2 . Since 
then other alloys have been discovered to possess the same attractive functionality. Although, 
the particular alloy, named Nitinol (Nickel, Titanium, Naval, Ordinance, Laboratory), is 
considered to be the SMA with the most potential for engineering applications because of its 
ductility at low temperature, high degree of shape-recovery capability (8% strain), large 
pseudoelastic hysteresis, corrosion and fatigue resistance, biomedical compatibility, and 
relatively high electrical resistance 3 . 

SMA’s are defined as a class of materials that have the ability to recover a given 
defonnation when heated above a certain temperature. There are two stable phases of the 
lattice structure for these metal alloys that may exist, tenned martensite and austenite. 
Austenite phase represents the high temperature stable state, and possesses the symmetric, 
cubic crystal lattice structure (See Figure 1.1) with a high elastic modulus. Martensite 
represents the low temperature stable state, and possesses the less symmetric lattice structure of 
a monoclinic crystal with a lower elastic modulus. Additionally, the martensite phase may 
exist as twinned or de-twinned (M + , M') states (See Figure 1.1). The interactions of thennal 
and applied mechanical loads to this material system are the basis for driving the phase 
transitions. The phase transitions between these two states are what provide the unique and 
attractive functionalities of SMA. When cooled under no applied mechanical load the material 
transforms from the austenite phase to the self-accommodating martensitic twinned phase. 
There is no observable shape change associated with this transfonnation. The boundaries 
between the layers of twinned martensite are easily deformed under applied load; therefore 
defonnation can easily be achieved in the martensitic state. Figure 1.2 shows the 
transformation behavior under heating/cooling and loading/unloading of an SMA element. 
Once a tensile load is applied in the martensitic state, the self-accommodated martensite- 
twinned configuration will transform, at a critical load, to M + or M' variants depending on the 
direction of the applied load. When the uniaxial load is removed the M + or M‘ phase will 
remain, leaving the given defonnation to the SMA element. Upon heating of the martensitic 
phase the “reverse transformation” occurs in which the lattice structure returns to austenite and 
recovers any defonnation. As previously stated this is tenned the shape memory effect. The 
SME always refers to the free recovery of the defonnation. It is important to mention that 
during the attempt to transfonn to austenite while constrained, high recovery forces may be 
generated. This makes SMA’s very attractive as actuators in mechanical systems. 

Associated with the martensitic and reverse transformations are four characteristic 
temperatures. Martensite start (M s ) is the temperature at which the SMA begins transfonning 



from austenite to martensite; martensite finish (Mf) is the temperature at which the SMA 
becomes fully martensitic and the transformation is complete; austenite start ( A s ) represents the 
temperature at which the “reverse transformation” begins; and austenite finish ( Aj ) represents 
the temperature at which the SMA is fully in austenitic phase. The difference in the phase 
transformation temperatures results in an associated hysteresis curve, which can be seen in 
Figure 1.3. 

SMA’s can also exhibit transfonnation pseudoelastic behavior, which is associated with 
the recovery of the strain upon unloading. If the material is above Af temperature and is loaded 
to a critical stress then the martensitic transfonnation can be induced, which is termed stress 
induced martensite (SIM). SIM produces a large defonnation in the material, which is fully 
recovered upon unloading when the material is above Af temperature (Superelasticity). 
Pseudoelastic behavior of SMA is attractive for use as an energy dissipation mechanism, and 
can provide hysteretic damping to structures. The pseudoelastic behavior of SMA is 
commonly exploited in medical and dental applications. 

SIM can also occur when the material is at lower temperatures than A/temperature, but 
only partial strain recovery occurs upon unloading at temperatures between A s < T <A/. In 
Figure 1.4 the stress-strain diagrams for various temperatures demonstrate the different effects 
of SMA phase transitions. When the material is at a temperature below Mf and is in the 
twinned martensitic phase its yield stress to induce fully M + de-twinned state, and thus a large 
defonnation, is quite low. Once the material is fully de-twinned and defonned to M + , 
continued loading would produce elastic defonnation until slip occurs. Once the material is at 
a higher temperature in the austenitic phase, the yield stress value is much higher than the 
martensitic phase’s. Figure 1.5 shows the conesponding load-defonnation curves for the 
material at various temperatures representing the quasi-plastic behavior (low temperature) and 
the pseudoelastic behavior (high temperature). 

As previously mentioned SMA’s are capable of producing large forces when 
recovering strain. The SME or reverse transfonnations described above refer to a free recovery 
(under zero load). If the strain recovery occuned under load the SMA performs work and is 
tenned a restrained recovery. If the SMA is prevented from recovering strain once heated 
above the transformation temperature that is tenned constrained recovery. Recovery stresses 
up to 700 MPa may be generated during constrained recovery. In this work the constrained 
recovery of SMA is used to induce large tensile forces in pre-strained SMA actuator ribbons. It 
is known that the recovery force generated during constrained recovery is a function of pre- 
strain, temperature, and thermal cycle. Cross reported the recovery stress vs. temperature for 
various pre-strain amounts, as seen in Figure 1.6 for Nitinol .1 in diameter rod 7 . At the A s 
temperature the recovery stress increases with temperature and then levels off near Af 
temperature. As pre-strain increases the recovery stress increases, but will never exceed the 
austenitic yield stress. Material characterization of Nitinol has been performed experimentally 
at NASA LaRC and Lockheed Martin Astronautics (LMA) to quantify the recovery stress vs. 
temperature, among the other temperature dependent properties, of the Nitinol alloy in ribbon 
form used for this study 8 . Figure 1.7 shows the recovery stress vs. temperature and thermal 
cycles at 4% pre-strain. It can be seen that the recovery stress/actuation degrades with 
thermal cycle to near 50 cycles and then levels off. 

Factors such as Ni content, aging, thermo-mechanical treatment and addition of 
alloying elements, which affect the structure, are important for controlling the memory 
behavior 9 . The transfonnation temperatures for SMA (M s , Mf, A s , Af) can be altered by 


2 



changing the composition of the alloy. It was reported by Beuhler and Wang that the A s 
temperature for Nitinol can be varied from -50 °C (-122 °F) to 166 °C (330 °F) by changing 
the Ni content from 54.4% to 56.5% by weight as shown in Figure 1.8 10 . When a nearly 
equiatomic Nitinol (Ni-Ti) alloy is fully annealed the SME consists of transfonnation directly 
from the austenite phase to the martensite phase. However, when the equiatomic Nitinol is 
partially annealed, the austenitic phase transforms to a R-phase, which consist of a 
rhombohedral lattice structure, before transforming to the martensitic monoclinic structure. 
This R-phase structure is an intennediate martensitic transfonnation, which is reported to 
strengthen the austenitic lattice structure and improve the recoverability of the SME 11 ’ 12 . This 
additional transfonnation, being a thennoelastic, martensitic transfonnation itself, also 
possesses a SME and pseudoelasticity but produces a much smaller strain recovery, and 
hysteresis. This is an attractive property for sensor applications. Due to the improved shape 
memory effect and strengthening of the material, partial annealing is desirable for most 
applications. The presence of the R-phase does not affect the perfonnance of the alloy for most 
applications, but will show up in material characterization measurements. 

Because there are so many aspects that may affect the transformation behavior of 
SMA’s it is important to have accurate material characterization measurements for experiment 
and modeling agreement. The transformation temperatures are typically measured by 
differential scanning calorimetry (DSC). This measurement technique results in plots of heat 
flow vs. temperature, and shows endothermic peaks corresponding to the austenitic phase 
transformations, and exothermic peaks corresponding to the martensitic twinned transfonnation 
and R-phase transfonnation. Thus, from the DSC measurements one can observe the values of 
the transfonnation temperatures, and the state of the SMA. Then mechanical properties of the 
SMA are measured by thermomechanical testing to quantify properties such as recovery stress, 
elastic modulus, and coefficients of thermal expansion vs. temperature to accurately model 
experimental behavior of SMA. 

1.2 Motivation/ Applications 

The unique functionality of SMA’s described above make them a very attractive 
material for a variety of applications. Since their discovery, early applications involved 
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coupling devices, fasteners, and electrical connectors . The most well documented large- 
scale application of SMA’s in 1971 was as a coupling to connect titanium tubing in the 
Grumman F-14 aircraft 9 . More recent applications for SMA’s have been focused in the 
medical field, and smart material research. “Smart Structures” have been defined by C.A. 
Rogers as material systems with intelligence and life features integrated in the microstructure 
of the material system to reduce mass and energy and produce adaptive functionality 14 . 
Materials commonly classified as “smart materials” include piezoelectric materials, 
magnetostrictive materials, electrorheolgical fluids, magnetorheological fluids, and shape 
memory alloys. All of these materials possess the ability to change their inherent properties as 
a result of a sensing mechanism. SMA’s have the ability to change stiffness, shape, natural 
frequency, damping and other mechanical characteristics in response to a change in 
temperature. The recent interest in smart structures has shown that SMA’s have great potential 
for vibration, structural acoustic, and structural shape control 15 . Smart materials such as 
SMA’s can be used, as actuator/sensing devices in structures to accomplish desired adaptive 
properties. One such type of structure was discovered when Rogers and Robertshaw 16 
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introduced the idea of embedding SMA actuators in a composite laminate. These structures 
became known as Shape Memory Alloy Hybrid Composites (SMAHC). 

The effects of noise and vibration are present in almost every aspect of every day life. 
It has become common practice to find ways to reduce noise radiation and vibration of 
structures. Common noise and vibration reduction techniques involve increasing the mass, 
stiffness, or damping of a structure. This is typically done by passive treatments such as 
constrained layer damping, acoustic absorption foam, tuned vibration absorbers, visco-elastic 
vibration isolators, etc. 17 . The concept of employing smart structures that can adapt to the 
surrounding enviromnent, to provide such previously mentioned functionality would be 
extremely attractive in any field of applications where weight-penalty is of concern. Once such 
application is in the field of high-performance aerospace vehicles where the excitation levels of 
their structures are high due to engine noise and turbulent boundary layer fluctuating pressures. 
Typical treatments for these effects are often limited to relatively low temperatures and high 
frequencies, and cannot adapt to a changing environment. The elevated thermal environment 
that aerospace vehicles are exposed to could activate SMAHC structures adaptive properties, 
rather than employ auxiliary hydraulic/electronic systems to facilitate actuation. Therefore 
smart structures such as SMAHC’ s are of great interest. 

There has been much research done in recent years on SMAHC’s. At NASA LaRC a 
great deal of experimental work has been perfonned on SMAHC structures with embedded 
Nitinol ribbon in which a new thenno-elastic model to predict the thennal buckling and post- 
buckling response of the structures was formulated 18 . Similar theoretical and experimental 
work has been done by Baz, who used Nitinol SMA wires for controlling the natural 
frequencies of composite beams 19 . The same has also been done for SMAHC plates 20 . A 
multitude of studies have been perfonned on the thennal buckling and post-buckling behavior 
of composite beams and plates reinforced by SMA components' ’ ’ ’ . Other work done in 
the area has included active vibration control of cantilevered beams with the actuators external 
to the structure' . There have also been analytical fonnulations developed to demonstrate the 
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vibration/structural acoustic control of SMAHC panels ' ’ . 

There has been much work done to develop theoretical models that capture the thermo- 
mechanical behavior of SMA, as well as characterization of SMA material properties ’ ’ ’ . 
There have also been many good comprehensive overviews of SMA characteristics and their 
applications published 6,9 . Although, as explained in the subsequent section, constitutive 
models that lend themselves easily to actual engineering property measurements are scarce. 

1.3 Objectives 

There has been previous investigation of SMAHC beam structures at NASA LaRC, 
which provided comparison of experimental results with a research code making use of a new 
thermomechanical model for SMAHC structures. The previous work showed qualitative 
validation of the model via the research code, but more definitive investigations are necessary. 
The goal of this work is to conduct experimental investigations on SMAHC beams to more 
accurately quantify the static buckling/post-buckling and dynamic response characteristics of 
the specimens. Numerical results will be generated in the commercial finite element code 
ABAQUS for comparison to experimental results, in contrast to the research code previously 
used. Fabrication of SMAHC beam specimens was previously perfonned at NASA LaRC, in 
which pre-strained (4%) SMA (Nitinol) ribbons were embedded within layers of glass-epoxy 
matrix material. The glass-epoxy matrix material was chosen in order to allow electrical 
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isolation of the SMA actuators, and flaw detection due to the material’s translucency. Once the 
SMA actuator ribbons are placed in the desired layers, the laminate is cured, and the SMA and 
glass-epoxy are integrated into one hybrid composite structure. Activation of the SMA 
actuators by resistive heating will serve as an adaptive stiffening mechanism to the structure by 
generating large recovery forces in the embedded pre-strained SMA ribbons. Accurate 
measurement techniques will be used to capture the adaptive stiffening effect of the SMA 
actuators. Numerical modeling will be performed in order to show the correlation between the 
measured and predicted responses. 

Results from measurements of the static response of the SMAHC beam specimens, 
under thermal load by resistive heating, will be presented. This requires the application of an 
electrical current through the embedded SMA ribbon. The A s temperature (40 °C (104 °F)) for 
the embedded SMA actuators is above ambient temperature thus once the temperature initially 
increases the SMAHC beam will be subjected to in plane compressive forces due to thennal 
expansion thus inducing buckling and post-buckling phenomena. The in plane compressive 
forces responsible for the post-buckled deflection will be overcome by the adaptive stiffening 
of the actuators once the temperature is high enough for activation of the SME. This highly 
nonlinear response behavior is quantified by a laser displacement transducer, Projection Moire’ 
Interferometry (PMI), and infrared thennal imaging measurements. These measurement 
techniques also accurately detennine initial geometric imperfections of the SMAHC beam 
specimens to help obtain accurate numerical models of the experimental response. 

Additionally, the steady-state dynamic response of the SMAHC beams, under random 
inertial loads is measured. Once the dynamic response time histories are captured for the beam 
specimens this data is post-processed into FFT, FRF, PSD, and coherence functions, which are 
used for damping estimates, and displacement RMS values of the SMAHC ’s vibration 
response. The results demonstrate the adaptive stiffening characteristics of the SMA actuators, 
and provide accurate data to be used in the numerical simulation of the experiments. 

A new thenno-elastic model was recently developed, that makes use of engineering 
property measurement, by Turner . This effective coefficient of thermal expansion model 
(ECTEM) is relatively simplified, and easily integrated into commercial structural analysis 
codes such as ABAQUS. The other aforementioned constitutive models are based on the 
micro-mechanical behavior of the SMA and typically calculate the stress in SMA material 
based on three components; an elastic component, transformational component, and a thermal 
component, that are all dependent on a martensitic volume fraction. The calculations of the 
volume fractions differ from model to model but are usually functions of stress and 
temperature. In the ECTEM the stress in an SMA material is based purely on an elastic 
component, and an effective thermal strain component. This effective thermal strain tenn 
represents both the thennal and transfonnational components of other models. As a 
consequence, the ECTEM model is limited to thennally activated transfonnations, but is 
particularly attractive for SMAHC because it only requires the experimental measurement of 
fundamental engineering properties. Therefore it can easily be used in commercial structural 
analysis codes that have classical nonlinear thenno-elasticity capabilities. In the case of 
SMAHC structures, the micro-mechanical models are not efficient at modeling many spatially 
distributed SMA actuators. The ECTEM allows this to be done easily, by implementing the 
composite properties from measurements of the material properties of each individual material 
making up the composite. The fonnulation of the ECTEM model will be explained in detail in 
Chapter 5 in order to show how a model making use of thermo-elastic properties can be 
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implemented easily into a finite element formulation for efficient structural analysis. Also, the 
reader is referred to the more detailed formulation. Turner 32 . 
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Chapter 2 


EXPERIMENTAL MEASUREMENT OF THE STATIC THERMOELASTIC 
RESPONSE OF THE SMAHC BEAM 


In this chapter a background is given on the fabrication of the SMAHC beam 
specimens, the buckling/postbuckling behavior, and geometric imperfections. Then the 
experimental set-up and instrumentation functionality are described. Then a brief overview of 
the Projection Moire’ Interferrometry (PMI) measurement technique is given. Finally, the 
experimental results of the laser displacement measurement and PMI measurements are 
presented and compared. 

2.1 Background 

The project consists of static/dynamic response experiments on two SMAHC beam 
specimens, which were previously fabricated. The fabrication procedure will briefly be 
described. For a more detailed explanation of the fabrication process the reader is referred to 
previous work at NASA LaRC 32 . The beam specimens were fabricated from glass-epoxy pre- 
preg and Nitinol ribbon materials with a ply orientation sequence of (45/ 0/-4 5/ 9 0)2 S - The glass 
epoxy pre preg has a thickness of 1.27e-4 m (.005 inches) and the Nitinol ribbon had cross 
section dimensions of .0023x1. 5e-4 m (,09x.006 inches). Therefore, 5 pre-strained (4%) SMA 
ribbons were placed side by side in the middle of the 0° glass-epoxy layer to represent the 0.45- 
inch actuator width called for by design. The schematic of the SMAHC beam specimen can be 
seen in Figure 2.1 with dimensions of 0.5588x0.0254x0.0019 m (22x1x0.078 in). The SMA 
actuator ribbons had a length of 0.6604 m (26 in), therefore there is an excess of SMA actuator 
material extending from the boundaries of the beam used as electrical leads for resistive 
heating. 

The SMAHC fabrication procedure consisted of pre-straining the SMA actuator 
ribbons (4%), stacking of each glass-epoxy ply in the corresponding orientation while placing 
the actuators in the 0.45 inch width strips cut from the middle of 0° plys, curing the laminate 
while restraining the actuators, and machining of the two beam specimens from the 
consolidated SMAHC laminate structure. The curing cycle consists of placing the laminate 
lay-up in vacuum drawn to an internal pressure of .0475 atm, and then subjected to an 
autoclave cure cycle in which a max temperature of 176.7 °C (350 °F) is reached. Upon curing 
of the lay-up, the laminate constituents are integrated into one hybrid composite structure. The 
whole hybrid composite structure panel and the beam specimens after machining can be seen in 
Figure 2.2. 

Typically when a beam is subjected to compressive loads in the axial direction, the 
beam loses stability. The Euler Load is the critical load at which the beam loses stability. 
Theoretically a lateral load is applied to cause deflection to a new position from which the 
beam will not return when the lateral load is removed. In reality this deflection occurs not as a 
result of an applied lateral load, but to geometric imperfections in the beam. The SMAHC 
beams were installed into a beam fixture with clamped boundary conditions. The first 
experiments will involve applying only a thennal load to the SMAHC beam by electrically 
activating the embedded SMA actuators. The heating of the constituent materials with 
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clamped boundary conditions will generate a compressive load within the SMAHC beam. 
Once this compressive load reaches the critical value (Euler buckling load) the beam buckles 
undergoing a large deflection. In classical theory the stiffness of the beam goes to zero and 
deflection goes to infinity at the critical load (P cr ). Mathematically the solution path of the 
beam “bifurcates” to a new equilibrium position. 

As previously mentioned the laminate consists of geometric imperfections, which 
dictate the direction of the buckling phenomena. As temperature of the post-buckled beam 
continues to increase, the embedded SMA actuators will be activated and generate in-plane 
tensile recovery forces that will overcome the compressive forces and render the post-buckled 
state of the beam back to its flat position. Therefore it is important to study these static 
responses of the SMAHC beams because the buckling/post-buckling states are part of the 
dynamic response as well. During the static response experiments accurate measurement 
techniques using Projection Moire’ Interferometry, laser displacement transducer, and infrared 
thermal imaging will be conducted to accurately capture the static response characteristics. 

2.2 Experimental Set-up 

In these static response experiments a thennal cycle, which consisted of increasing the 
temperature of the SMA actuators from 23°C (74 °F) to 121 °C (250 °F), was conducted 10 
times for 2 separate beam specimens for repeatability investigation. The SMAHC beam was 
installed into the beam fixture as shown in Figure 2.3 with clamped boundary conditions. Each 
SMAHC beam has an unsupported length of .4572m (18 in). The aluminum beam fixture 
consists of a mechanical grip, which applies fixed boundary conditions on each end, and 
electrical grips that apply the electrical current to the SMA actuators embedded within the 
beam. Figure 2.4a shows the mechanical and electrical grips from the front of the beam 
fixture. The positive power supply lead can be seen connected to the electrical grip block. 
This entire block is electrically activated but is isolated from the mechanical gripping block by 
a layer of Garolite G-10 insulation. The SMAHC beam is also thermally isolated from the 
beam fixture by layers of Garolite G-10 installed between the mechanical grip surfaces. The 
Garolite G-10 is a glass fabric-epoxy insulator, which is a relatively hard/stiff material that 
should eliminate any compression under in-plane loads caused by the SMA recovery stress 
during testing. In similar tests a fibrous ceramic insulator was in service, which showed signs 
of compression during testing, thus relieving some of the recovery stress in the activated SMA 
actuators. Type-T thennocouples were attached to the beam specimen on the center of the 
beam 1 inch from the left mechanical grip, to both of the mechanical grips, and one is 
measuring the surrounding ambient air temperature. These thennocouples can be seen attached 
to the left mechanical grip and the beam specimen in Figure 2.2b, which shows the backside of 
the beam fixture. Thennocouple measurements are subject to large error due to contact 
resistance, and “fin” effect, in which the attached thennocouple lead conducts heat away from 
the contact and the heat is lost into the sunounding air due to convection. Therefore, the 
thermocouple measurements are used for monitoring purposes only. An infrared thermal 
imaging camera (FLIR Thermacam SC2000) was used for control measurements of the beam 
specimen temperature. The thennal camera eliminates the measurement uncertainty of contact 
resistance and “fin” effect, but does require calibration to the target emissivity. The average 
total hemispherical emissivity of the beam specimens was measured to be 0.8 and it was 
assured that the beam radiates as a gray-body. 
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A Hewlett Packard DC Power Supply, model #6652A, is connected to the two 
electrical grips of the beam fixture, and is controlled via a LabView virtual instrument (VI). 
The LabView VI controller was previously configured and designed at NASA LaRC. The VI 
controller makes sensor measurements through the infrared thennal imaging camera, the type- 
T thermocouples, and controls the power supply current output through general-purpose 
interface bus (GPIB) commands. Optimal gains are set to use proportional-integral-derivative 
(PID) control to supply current to raise the temperature of the SMAHC beam centerline to 
desired set points where data is recorded. The LabView VI system is capturing data from a 
laser displacement transducer (Keyence model LK-081). The laser displacement transducer 
was attached to an optical positioning traverse, which was used to precisely align the laser to 
the mid-span point of the SMAHC beam specimen, and calibrate the initial laser 
transducer/beam surface distance to zero. This was done before every run. Therefore as the VI 
controller runs through a thennal cycle, the beam temperature is sensed by the thennal camera, 
cunent is applied to raise the temperature from 23 °C (74 °F) (ambient) to 121 °C (250°F) 
capturing set points in increments of every 5.55°C (10 °F), and the mid span displacement 
values are being recorded continuously approximately every 4.2 seconds during the thennal 
cycle except for beam 3 run 2 and beam 3 run 3, in which laser data was captured at set-point 
only. One entire cycle took approximately 480 seconds. As an added feature the VI controller 
was designed to automatically trigger the PMI imaging system to capture images at every set 
point during a thennal cycle. Therefore, the laser displacement transducer is capturing the 
mid-span displacement measurements, and the PMI system is capturing the full-field 
displacement of the SMAHC beam specimen, throughout the thermal cycle. Figure 2.5 is a 
schematic showing an overview of the components and functionality of the experimental set- 
up. Figure 2.6 is an image of the entire experimental set-up and all measurement components 
can be seen. The PMI measurement of full-field displacement will be explained in more detail 
in the subsequent section. 

The results of the full-field displacement measurement give great insight into the 
buckling behavior of the beam, and initial imperfection measurements can be captured. Also 
the PMI displacement measurements can be compared to the laser displacement transducer 
measurements of the SMAHC beam mid-span to reduce uncertainty. 

2.3 Projection Moire’ Interferometry 

Projection Moire’ Interferometry is an optical surface defonnation measurement 
technique used to measure the full-field deformation of a structure. PMI is an optically simple, 
non-contacting measurement technique used since the 1970’s for surface topology and shape 
characterization 33,34 . A projection system with an installed Ronchi ruling is used to project a 
grid of equally spaced, parallel lines onto the structures surface. A CCD camera is positioned 
to view the structure at some angle 9 inclined from the projector optical axis, which is 
perpendicular to the structure surface. Images are captured of the projected grid lines during 
the initial unloaded state and during loading of the structure. As previously mentioned for this 
experiment, images were captured every 5.55°C (10°F) from 23-121°C (74-250°F). At each 
set-point the PMI system is triggered, and it takes 10 images and averages them. Image 
processing routines are used to remove camera perspective distortion and interfere the captured 
images with a numerically generated reference grid. This results in interferograms containing 
moire’ fringes. The fundamentals of PMI measurement and its applications have been studied 
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and quantified at NASA LaRC by G.A. Fleming for use in wind tunnels . Figure 2.7 shows a 
schematic of the PMI imaging system. 

Imagine that the grid lines are projected onto a flat surface located at coordinate z 0 in 
the schematic. The CCD camera captures this reference image. Then when the plate is 
deformed or translates in the z direction the spatial locations of the gridlines change, thus 
producing a deformed image. This change can be represented mathematically by distorting the 
reference image grid by this phase function: 


<p(Uj) 
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( 1 ) 


where <p(i,j) is the spatial phase at a pixel in the CCD camera image plane corresponding to 
point (x,y) on the un-deformed object surface, z(x,y) is the out-of-plane deformation at that 
point, and d is the distance between grid lines. It is known that taking the absolute difference 
of the reference grid and the deformed grid produces a interferogram. The interferogram will 
contain moire’ fringes occurring when cp{i, j) = Inn , thus: 
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Therefore, the moire’ fringes in each interferogram represent the contours of the defonnation in 
the z direction. The fringe contour interval is then defined as: 
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d 
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This is the basic idea of how the PMI measurement theory is implemented. Further 
techniques of image processing such as phase shifting, phase unwrapping, choosing anchor 
point, and grid-line nonlinearity issues are addressed in more detail in the following 
references 35,36 ’ 37 . 


2.4 Static Thermoelastic Response Experimental Results 

Laser Displacement Transducer Results 

The mid-span deflection measurements of the SMAHC beam2 and beam3, by the laser 
displacement transducer, are presented in Figure 2.8, 2.9, respectively. The figures show the 
mid-span deflection as a function of average temperature of the centerline along the length of 
the beam. Data from thermal cycles 1-9 on SMAHC beam specimen #2, and thennal cycles 1- 
10 on SMAHC beam specimen #3 are presented. The adaptive stiffening mechanism is evident 
in all runs for both beams. Recall that the austenite start temperature for the Nitinol ribbon is 
greater than the ambient temperature (A s > T 0 ). Therefore, buckling of the clamped beam will 
occur immediately when the temperature of the beam begins to change. As the temperature of 
the beam is increased the midspan of the beam deflects until a critical temperature at which the 
SMA recovery forces overcome the in-plane compressive forces due to thennal expansion and 
the beam is rendered flat. The experimental data shows a relatively large variation between 
maximum buckling deflection. This is a highly nonlinear phenomenon and dependent upon 
many parameters such as initial deflection, initial stress, geometric imperfections, and 
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boundary conditions. Each thermal cycle test and possible reasons for the behavior will be 
described in the next paragraph. 

In Figure 2.8 for beam 2, runs 1-9 are the first 9 thennal cycles on this beam specimen. 
In Figure 2.9 for beam 3, runs 1-10, correspond to thermal cycles 5-14. Actuator degradation 
with thermal cycle reduces the maximum tensile recovery force in the beam at higher 
temperatures. Figure 4.4 shows the measured recovery stress as a function of temperature for 
numerous thennal cycles of the pre-strained Nitinol ribbon used in this project. Notice, even at 
high thermal cycles the temperature range that corresponds to where the post-buckled 
deflection is recovered (~32 °C - 65.5 °C) shows no significant degradation with thermal cycle. 
Therefore thermal cycle degradation is not a factor in the observed variation of maximum 
deflection. In run 1 for beam 2, and runs 1, 2, and 10 for beam 3 a relatively large max 
deflection occurs. All these runs occurred as the first thennal cycle on the SMAHC beam after 
a relatively long period of storage time, and could be attributed to moisture absorption in the 
beam. It is common for glass-epoxy composite materials to exhibit moisture absorption after 
relatively short periods of time. The absorption of moisture in the composite beam laminate 
could possibly cause swelling of the layers of the structure. The presence of moisture in the 
composite structure could enhance the thermal expansion effect at the initial temperature 
increase. As temperature continues to increase the maximum buckling deflection is increased 
and the point of deflection recovery is pro-longed to a higher temperature. It is likely that at 
some temperature the moisture is driven out of the structure, thus causing a rapid recovery of 
the deflection. This could be the explanation for the rapid deflection recovery observed in the 
large deflection cases. Beam 3 Run 1 was the first thennal cycle on the beam in approximately 
one year. Beam 3 Run 2 occurred 9 days after the first run. Beam 3 Run 10 occurred 14 days 
after the previous run. Beam 2 Run 1 was the first thermal cycle on the beam, which was 
fabricated approximately one year before the experiment. So, it seems probable that the 
buckling phenomenon and deflection behavior are very sensitive to hygothermal effects, and 
indicates that the numerical modeling to match the buckling phenomenon will require precise 
knowledge of all parameters that can affect the environmental state of the beam. 

Another contributing factor to the variation could also be a result of pre-tension on the 
SMA actuator leads once installed. It is important to have no stress present in the installed 
SMAHC beam initially, but it is also necessary to have the beam installed “taught” so the beam 
is not initially buckled. The experimental procedure involved pulling the leads “taught” just 
before clamping at the boundaries. It is a possibility that some pre-stress may be present in the 
SMA actuator ribbons initially before the run. Recall that the transfonnation temperatures are 
a function of stress, and typically vary linearly with applied stress. The buckling deflection is 
very sensitive to the recovery force of the SMA, which would be effected by the shift in the A s 
transformation temperature. During the buckling and post-buckled state of the beam, there are 
two mechanisms competing. The compressive forces being generated due to thermal 
expansion, and the recovery stress that is generated at A s temperature of the embedded SMA 
actuators. A small change in either mechanism could result in a large change in the deflection 
behavior. 

It should also be mentioned that another possible contribution to the pre-stress variation 
could perhaps be due to small deformation of the Aluminum beam fixture in which the 
SMAHC beam is clamped within. During the experiment procedure thermal cycles were 
perfonned in succession, once the infrared thermal imaging camera measured the beam’s 
centerline cooling from elevated temperatures back to ambient (~74 °F). In cases, where the 
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thermal cycles were performed in quick succession, the Aluminum fixture may have been 
somewhat soaked with thennal energy, which it looses rather slowly relative to the SMAHC 
beam specimen. As a result, when the SMAHC beam was re-gripped during these cyclic runs 
the fixture could be continuing to lose heat and contract. Thus, the beam would be under 
compressive stress at the beginning of a thennal cycle, which was perfonned in quick 
succession from a previous one. 

Another observation from the laser displacement data is the variation of the initial 
imperfection of the mid-span. Beam 2 seems to exhibit a larger variation in initial deflection 
values between runs, than does beam 3. This could be attributed to the fact that beam 2 had 
never been subjected to a thennal cycle before, thus the internal stress state probably has a 
large variation in the first couple of thennal cycles. After curing of the laminate, residual 
stresses are present in the composite matrix. Upon the first thennal cycle the largest recovery 
forces are seen within the beam, causing a large chance of internal slipping, and settling of the 
composites components during early runs. Thus the initial deflection of the mid-span, which is 
due to geometric imperfections of the composite, and also induced by the boundary condition 
imperfections of the fixture, could have larger variations if the internal state of the beam is 
changing. Another interesting observation is that during all runs for both SMAHC specimens 
the buckling phenomenon always occurred in the same direction towards the laser 
displacement transducer. For Beam 2 run 9 the SMAHC specimen was flipped over in the 
beam fixture, and still exhibited buckling towards the laser. This suggests that boundary 
conditions of the beam fixture are contributing to the initial shape, and buckling direction of the 
beams. 

Following the static experiments in which the laser displacement data and PMI data 
was captured, another experiment was run to capture the geometric imperfections of the beam 
and capture the full-field initial shape of the SMAHC beam installed in the fixture to confirm 
the PMI results. The experiments uncovered some of the initial conditions contributing to the 
initial state of the beam. 

The experimental configuration consisted of a Klinger Scientific motorized linear 
traverse and controlled via the Labview VI. The translation stage was mounted parallel to the 
SMAHC beam fixture with installed beam as seen in Figure 2.10. The laser displacement 
transducer was mounted to the translation stage, and traversed via a programmable stepper 
motor the length of the beam and recorded displacement measurements in .10 inch increments 
along the 18 inch length of the beam. This data shows the initial imperfection of the beam, 
when mounted. The laser is calibrated to zero at the left end of the beam from a front view, 
and all subsequent measurements during a run are relative to that first set point. Multiple cases 
were studied, and these experiments uncovered some interesting results. 

The initial imperfection data was captured on the beams in one direction and then 
flipped to the other side in the beam fixture. These results in Figure 2.11 show that the 
SMAHC beams take on the same shape even once flipped. This figure shows the out of plane 
deflection of the beam at ambient temperature as a function of the length of the beam. These 
measurements represent the initial imperfection profile of the beam when installed in the 
fixture. Notice that the blue colored line represents the out of planed deflection measurement 
for beam 2. Then beam 2 was flipped over and installed again in the fixture, which it would be 
expected to see the opposite deflection profile. But the magenta line results for Beam 2 flipped 
produces the same deflection profile. This confirms that the boundary conditions are 
contributing to the initial shape of the beam. Approximately the same deflection profile results 
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for beam 3 specimen as well. The laser data is also compared to one of the data sets for the 
PMI absolute imperfection measurement. As will be explained in the next section the absolute 
PMI measurement is the first image of the beam captured at ambient differenced from a 
reference image of a perfectly flat calibration plate. It is evident in Figure 2.1 1 that the laser 
displacement initial imperfection results are in excellent quantitative agreement with the black 
dashed line representing the PMI absolute measurement at ambient temperature. This is true 
for all runs of the PMI absolute measurements for both beams at ambient and confirms that the 
PMI absolute measurements are correct. 

An interesting observation is that the beam is skew in all cases from one grip to the 
other in Figure 2.1 1. At 18 inches represents the left boundary grip therefore it is expected that 
deflection measurement should be zero relative to the other right boundary grip This was 
noticed in the PMI results and first attributed to the possibility of the calibration plate not being 
perfectly aligned with the beam fixture. But, the laser displacement measurements confirm 
that there is an offset in the actual fixture from one grip to the other while a beam is installed. 
A measurement of the distributed thickness of the Garolite G-10 insulation blocks that the 
beam is clamped within show that they are not symmetric. In fact the difference in thickness of 
the inside edges of the two back G-10 insulation blocks that the beam is clamped against have 
a difference of approximately .020 inches. This difference closely corresponds to the skewness 
observed in the PMI, and laser displacement measurements of the initial imperfection. 
Subsequent measurements using the laser displacement transducer traverse with multiple G-10 
block arrangements were perfonned to show that the effect the boundary conditions have on 
the initial shape of the beam. The results in Figure 2.12 show a totally different initial shape of 
the beam for different orientations. The corresponding configurations are best described by the 
schematic showing the G-10 block arrangement and each G-10 thickness at the edges in Figure 
2.13. 

Thus, based on the above observations the laser displacement transducer translation 
experiments prove that the boundary conditions are the primary factor to the initial shape of the 
beam, and confirm the PMI absolute measurements at ambient. The offset of the beam 
between boundaries is a physical entity of the experimental results, but is just a rotation of the 
beam. All numerical modeling in ABAQUS will be relative to a flat plane, so for comparison 
purposes all skewness present in the initial geometric imperfection experimental data has been 
subtracted before input into the ABAQUS input file. 

Projection Moire ’ Interferometry Results 

The PMI measurements captured the full-field displacement of the beam and are 
represented as pixels. The full-field PMI images showed no variation in displacement in the 
vertical width direction, thus the centerline data points along the length of the beam are 
representative of the full-field image of the beam. The PMI data presented for this experiment 
shows the out of plane displacement of the centerline of the SMAHC beam for various SMA 
actuator temperatures. The x-axis represents the beam length. PMI data points were only 
captured for 0 - 17.83446 inches along the length of the beam due to the angle of the PMI 
projector creating a shadow at the right boundary. PMI data is relatively noisy due to surface 
roughness, and processing techniques. To adjust for this, experimental data was fit to a 6 th 
order polynomial of the fonn: Y = A + Bx + Cx 2 + Dx* + Ex 4 + Fx 5 +Gx 6 , in which data 
points from -1.006 to 0 and 18 tol9.012 inches were manually zeroed to produce accurate line 
fits representing the defonnation profile of the beam. The zero padded data sections represent 
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the x-locations in the clamped boundaries and it is a fair assumption that displacement values 
are zero. Thus, the line fit was used to extrapolate data for the small beam length section from 
17.83446 - 18 inches. The positive sense defonnation represents deflection away from the PMI 
imaging system. An example of the captured PMI data and corresponding line fit for a 
differential measurement at 48.8 °C (120 °F) is presented in Figure 2.15 for Beam 2 Run 4. 

There are two types of PMI measurements made for these experiments. As explained 
the PMI data can only represent deformation from an initial reference image. Differential 
measurements use the initial ambient temperature image as this reference image. This 
reference image is considered the zero state. Subsequent images taken as temperature 
increases are the deformations relative to the reference state. Therefore any initial imperfection 
of the SMAHC beam is not present in these measurements. Additionally, the skewness is not 
present in the differential measurements for the same reasons. 

The other type of PMI measurement is the absolute measurement. In reality it is known 
that the SMAHC beam specimen in the initial ambient temperature is not perfectly flat, due to 
imperfections in the composite laminate lay-up, initial stress state within the beam, and 
imperfect boundary conditions. The absolute PMI measurements use a reference image of a 
flat plate in the plane of the beam fixture (See Figure 2.14). The first ambient temperature 
image captured from the SMAHC beam is then differenced with the flat plate image. This 
results in a measurement representing the full-field initial imperfection of the SMAHC beam at 
ambient temperature. 

The Differential PMI measurements were post-processed by interfering the reference 
image at ambient with each elevated temperature image, which gives the difference. The fact 
that the reference image was skewed, as well as all elevated temperature images it is 
differenced with, results in no offset present in the differential measurements. Therefore the 
confirmed initial imperfection, with skewness subtracted off, can be added to each differential 
measurement to give the absolute shape of the beam for all temperatures in order to show 
comparison with the numerical simulations. Each run’s absolute measurement, giving the 
initial geometric imperfection, was added to each differential measurement. Therefore, there 
now is a representative absolute measurement for every temperature in 5.55°C (10°F) 
increments between ambient and 121.1 1 °C (250 °F). The PMI results are shown for selected 
runs 1, 4, and 8 for beam 2, and beam 3 respectively in Figures 2.16-2.21. The plots show the 
full field out of plane deflection of the beam for a few individually selected temperatures. 
These presented results in the aforementioned figures show the same trends as all runs, which 
are included in Appendix A. Notice that in all runs the initial shape of the beam takes on the 
same S-shaped geometric imperfection at ambient temperature. This will be the initial 
geometric imperfection data input into ABAQUS for accurate modeling of the experiments. 
Also notice that the adaptive stiffening is also evident in all the PMI data and matches the laser 
displacement results for the midpoint of the beam. 

The midspans of the centerline PMI deflection profiles were extracted for each 
temperature and plotted with the laser displacement transducer midpoint deflection 
measurements in Figure 2.22, and Figure 2.23 for beam 2 and beam 3 respectively. It is 
evident that there is good agreement between the two measurements. Note that there are fewer 
data points making up the PMI data curve because there is only PMI midspan data at the 
temperature set points. Also, small differences in the data sets can be attributed to the fact that 
the PMI data system was triggered to take data once the temperature reached a certain set point. 
The PMI system is triggered after the set point temperature is reached and laser data is already 
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written. Thus, 10 averages are collected as the thermal controller starts raising the temperature 
to achieve the next set-point. A similar scenario is likely in Beam 3 Run 1 because this is the 
only run where the control voltage was operated manually. Therefore there is a greater chance 
that the set point temperature drifted from the desired temperature by the time the PMI was 
manually triggered. This can explain the large difference in this particular data set. There are 
also small differences between data points at ambient temperature. This can be attributed to 
the fact that they are two separate measuring devices that are not referenced to the same initial 
position. The midspan data point for the PMI is referenced to the flat calibration plate image in 
the plane of the beam, whereas the laser displacement transducer’s initial data point is 
referenced to the beam’s flat position at 121.1 °C (250 °F). 

Overall the PMI measurements prove to be extremely accurate, and in excellent 
agreement with the laser displacement transducer measurements. This validates the accuracy 
of the PMI measurement technique for this experiment. The PMI measurements give 
additional insight into the full-field defonnation of the SMAHC beam throughout a thennal 
cycle, and give the important initial imperfection data to be used for accuracy in numerical 
modeling. 
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Chapter 3 


EXPERIMENTAL MEASUREMENTS OF THE DYNAMIC RESPONSE OF THE 

SMAHC BEAM 


The dynamic response chapter is arranged as follows. First, a background is given on 
the dynamic response measurement of the SMAHC beams. Then the experimental set-up and 
instrumentation functionality are described. Then a brief overview of the post-processing of 
the captured data is presented. Finally, experimental results of the SMAHC beam specimen’s 
responses are presented and compared for a repeatability investigation. 

3.1 Background 

All the SMAHC beam vibration results in this study result from random base 
acceleration excitation. Resulting responses will be presented at various temperatures 
corresponding to the different stiffness states during the thermomechanical response. Using a 
random inertial excitation requires a statistical analysis to be used to process the response data. 
Randomness in vibration analysis can be thought of as the result of a series of experiments all 
perfonned in an identical fashion under identical circumstances, each of which produces a 
different response 4 ". Multiple time histories of the input and corresponding responses are 
sampled. From these multiple samples Fast Fourier Transforms can be done and averaged to 
produce dynamic response characteristics. By demonstrating control over the dynamic 
response, i.e. adaptive stiffening, and RMS value reduction, implications of weight-efficient 
vibration and acoustic radiation reduction methods are eluded too. Also analysis of the 
resulting frequency responses will result in beneficial modal parameter values to improve 
modeling attempts. The beam fixture was designed not to have any resonance in the excitation 
bandwidth used in the experiment. Transducers are placed on the fixture and the beam in order 
to obtain the relative response of the SMAHC beam, and allow for comparison to response 
spectra generated numerically by ABAQUS. During excitation a thennal load will be applied 
to the beam as in the previous static tests. The SMAHC beam will be brought to a desired 
temperature set point, in which dynamic data will be captured. One run will consist of 
capturing dynamic data at all the chosen set points through a thennal cycle from ambient to 
121.11 °C (250 °F). 

3.2 Experimental Set-up 

The aluminum beam fixture with installed SMAHC beam specimen is mounted to a 
large electromechanical shaker (Model T1000) configured horizontally and attached to a 
magnesium slip table, which oscillates over a bed of oil in a horizontal plane during operation, 
as seen in Figure 3.1. The excitation of the beam during operation is measured by 3 Endevco 
Model 2250 A- 10 accelerometers mounted in 3 different positions along the length of the 
SMAHC beam specimen. There are also 2 PCB model T352M92 accelerometers mounted to 
each end the beam fixture itself. Two of the PCB accelerometers are used for control of the 
shaker table, and the other two are used to subtract off the base motion from the absolute 
motion measured by the accelerometers mounted to the beam. The dynamic response 
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experiments of each SMAHC beam specimen consisted of 3 thermal cycles where the beam 
was heated from ambient to 121.11 °C (250 °F) in increments of 5.55 °C (10 °F). At each set 
point the dynamic acceleration data was captured. The experimental set up can be seen in 
Figure 3.1. The infrared thermal imaging camera is aligned at the end of the shaker table 
focused on the SMAHC beam, and used for control. A DC power supply voltage is applied to 
the beam via the Labview VI controller in the same manner as the static experiments described 
previously. The Labview VI controller is sensing measurements from thermocouples for 
monitoring, and the infrared thermal imaging camera for control. An addition to the Labview 
controller was made to hold once a set point temperature is reached, and then can be released 
when specified by the user. This allowed the dynamic data to be captured for that temperature, 
and then the set point is released and the beam’s temperature is elevated to the next set point. 

In Figure 3.2, there are Endevco (2250A-10) accelerometers mounted on the SMAHC 
beam in positions 1-3. Position 1 is 4.5 inches from the right mechanical grip from the back 
view, position 2 is at the midspan, and position 3 is 14.5 inches from the right mechanical grip. 
There are also 2 PCB (model T352M92) accelerometers mounted on the back of each of the 
mechanical gripping blocks. One pair is used for control of the shaker table, and the other pair 
is used to capture the fixture response. Time histories with a total length of 200 seconds were 
captured with a sampling rate of 1024 Hz to allow 50-frame averages with a bandwidth of 0- 
512 Hz and a frequency resolution of .25 Hz for every set point. The SMAHC beam was 
excited by random base acceleration over a 10-400 Hz bandwidth with a RMS value of .25 g. 
The excitation’s power spectral density (PSD) is shown in Figure 3.3. All three response 
transducers are capturing acceleration data, and one pair of the fixture transducers are being 
integrated to displacement data, while the other pair for control of the shaker is capturing 
acceleration data. All of the transducers are connected to a Hewlett Packard VXI mainframe 
(model E8403A) with a 16 channel E1432A module data acquisition card. The dynamic data 
is captured and processed using MTS IDEAS software. Upon capturing all of the time 
histories for each temperature set point, the data is post-processed into relative response 
quantities by subtracting off the base acceleration. Then response quantities are generated 
including frequency response functions (FRF) between the beam specimen’s response and the 
fixture’s input, power spectral densities (PSD) of the acceleration input and relative beam 
responses, and coherence functions, which indicate data quality in general. A schematic of 
the dynamic experimental set-up can be seen in Figure 3.4. 

Before the dynamic response experiment the beam specimen was weighed and installed 
into the beam fixture. Beam 2’s mass without mounted transducers was 75.97e-3 kg (.1675 
lbf) and 79.15e-3 kg (.1745 lbf) with mounted transducers. Beam 3’s mass without mounted 
transducers was 74.84e-3 kg (.1650 lbf), and 78.25e-3 kg (.1725 lbf) with transducers. Just 
prior to activating the excitation to the base, the mechanical and electrical grips are torqued to 
13.56 Nm (120 inlbf). Once the random broadband excitation begins the thermal controller is 
activated and holds at the designated set point. The dynamic data acquisition system then 
acquires the acceleration data for all accelerometers. Following three runs for each beam the 
acceleration time histories are post-processed, which included averaging 50 times. 

3.3 Dynamic Response Data Post-Processing Overview 

In general a transducer used to measure input and output during a vibration test usually 
contains noise. Additionally, perfonning a number of identical tests and averaging the results 
increase confidence in a measurement of random vibration. In this experiment it was desirable 
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to obtain displacement response data, because it is more intuitive. Direct hardware integration 
was attempted during testing by PCB model 480B10 Signal Conditioners. In this case 4- 
second time history samples were taken to construct a 200 second time history for each 
temperature set point in MTS IDEAS. The integrated time signals were swamped with a low 
frequency noise. The acceleration time signals were processed, and the FRF’s and PSD’s 
showed large low frequency peaks representing the noise. Therefore, upon integration of the 
acceleration data in the frequency domain the low frequency noise dominates the resulting 
displacement time history. Thus, it was determined that all transducers would capture the 
acceleration data, and integration to displacement would be performed offline. A MATLAB 
m-file was written to perfonn a Fast Fourier Transfonn (FFT) on the acceleration time history, 
integrate in the frequency domain to get displacement FFT values, “zero out” the low 
frequency noise below 10 Hz, and then perform an inverse FFT to convert to a displacement 
time history data file. Additionally this m-file computes the displacement PSD’s, and RMS 
values from the generated displacement time history (See attached m-file in Appendix B). 

The Fast Fourier Transform is defined as 

X{f,t)= ex p[- j2/fnAt] (4.1) 

n = 0 


assuming that there is a time history record x(t) that is sampled at N equally spaced points a 
with time increment At, where At has been selected to produce a sufficiently high cutoff 
frequency. Therefore, x n =x(nAt) where n = 0, 1, 2,...., N-l. The selection of the 
corresponding frequency values for the computation of X(f,t) is 

f k =~ = — k = 0,1,2,..., N-l (4.2) 

T NAt 


At these particular frequencies, the transformed values give the Fourier components defined as 




At 


n = 0 


. 2 nkn 


k = 0,1, 2,..., N-l 


(4.3) 


The Fourier components repeat after the index number k=N/2 since the Nyquist frequency 
occurs at this point. Thus, in the m-file the corresponding frequency vector generated for the 
Fourier components is folded back on itself in order to retain the phase infonnation to 
transfonn back to the time domain. 

Matlab generates the PSD’s using the specified time data, number of points per data 
segment, sampling rate, window vector, and number of points to overlap. The output of the 
PSD function in Matlab must then be nonnalized to the Nyquist frequency to obtain the proper 
units of (in)“/Hz. The RMS value is then calculated from the auto spectra (PSD) values. In 
Figure 3.4, a particular displacement PSD at ambient temperature generated by the m-file is 
plotted with its corresponding acceleration PSD to confirm that the m-file is correct. It is 
observed that the low frequency noise is filtered out, and the frequency infonnation is retained, 
providing excellent correlation with its corresponding acceleration PSD. 
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3.4 SMAHC Beam Dynamic Response Experimental Results 

The experimental behavior of the SMAHC beam specimen observed during dynamic 
testing will briefly be described. As described in chapter 3, the beam exhibits a thennal post- 
buckling deflection over a particular range of temperatures. Within this temperature range the 
beam undergoes interesting dynamic responses because the stiffness of the beam is changing 
rapidly. Recall from the post-buckling analysis in chapter 2 that the stiffness theoretically 
drops to zero at the point where the beam loses stability and buckles. When the beam is in the 
post-buckled state the stiffness increases with deflection. When the SMA actuators are 
activated and return the beam flat, the stiffness again decreases but rapidly increases with 
continued heating from this state. For accurate numerical predictions in this temperature range 
a nonlinear dynamic analysis is needed, and will prove difficult to get good correlation for the 
simple fact that the beam’s dynamic characteristics are changing rapidly. During excitation of 
the beam, its relative response was very noticeable at temperatures where the SMA actuators 
were not activated; where-as at high temperatures the relative response was small, and 
unnoticeable once the SMA ribbon is fully activated. Dynamic results will be presented for 
these temperature ranges. The SMAHC beam is returned to the flat configuration at 
approximately 65.55 °C (150 °F). 

In Table 3.1 and Table 3.2 the RMS displacement values are presented as a function of 
temperature for all 3 positions on beam 2, and beam 3 respectively. The results show a 
decrease in RMS displacement of the beam with increasing temperature, which is evidence of 
the adaptive stiffening of the beam. The modal parameter estimates for beam 2 and beam 3 
were obtained and are presented in Table 3.3 and Table 3.4 for each run. The modal damping 
ratio estimates are particularly useful as input into the numerical modeling of the experiments 
in ABAQUS. These modal damping ratio values are obtained in MTS IDEAS by perfonning a 
single degree of freedom curve fit (SDOF method). In this method each peak is analyzed by 
assuming that it is the frequency response of a SDOF system. This assumes that near each 
resonant peak, the FRF is dominated by that single mode. The damping estimates for mode 1 
and mode 3 show a decrease in damping from the non-activated state relative to the high 
temperature activated state of the SMA. Interestingly, at the higher temperatures near 121.1 °C 
(250 °F) the damping estimates begin to increase slightly. This may be attributed to the glass 
epoxy matrix material heating up and softening, and is not associated with the SMA actuator 
material. 

In Figures 3.6-3.11 the displacement PSD’s for selected temperatures of 23.3 °C (74 
°F), 71.1 °C (160 °F), 93.3 °C (200 °F), and 121.1 °C (250 °F) are presented for the various 
positions along the SMAHC beam specimens. The PSD’s for position 2 differ from the 
corresponding PSD’s for positions 1 and 3 in that the midspan of the beam does not contain the 
anti-resonant peaks between mode 1 and mode 3 resonance peaks. Also position 2 contains a 
larger displacement RMS value, which makes sense at the midspan of the beam. The 
experimental results show excellent repeatability between runs 1, 2, and 3, as well as between 
beam 2 and beam 3. Therefore the presented data sets in Figures 3.6- 3.11 are representative 
of all other captured dynamic data. The PSD’s for both beams at all positions demonstrate the 
stiffening effect of the SMA actuators as the modal frequencies of the response are increased 
with increasing temperatures above 60 °C (140 °F). Figure 3.12 best shows the variation of the 
modal frequencies with temperature. Notice that the fundamental frequency increases from 
30.1 Hz to 125.1 Hz from ambient temperature to 121.1 °C (250 °F). Also, mode 3’s frequency 
nearly is increased out of the bandwidth of excitation. 
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Thus the dynamic response experiments demonstrate the adaptive stiffening mechanism 
of the embedded SMA actuator ribbons. The dynamic data gives a good overview of the full 
dynamic behavior for multiple positions along the beam length. The fact that the RMS values 
decrease with temperature, and the beam is stiffened indicates the reduction in the dynamic 
response. This has implications of improving sonic fatigue, and noise transmission reduction. 
Finally, accurate modal parameter estimates are captured to include in the dynamic numerical 
modeling of the SMAHC beam in ABAQUS. 
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Chapter 4 


THERMOELASTIC CONSTITUTIVE MODELING FORMULATION 


4.1 Background 


32 

The ECTEM developed by Turner makes use of experimentally measured 
engineering properties. This makes the ECTEM easily implemented into a commercial code 
that is capable of nonlinear thennoelastic analysis of composite materials. This chapter will 
briefly summarize the development of this model, to give the understanding of how the 
engineering property measurement becomes the detennining parameter of the modeling. The 
model is based on a mechanics of materials approach and the reader is referred to numerous 


sources 


18,38,39,40 


4.2 Thermoelastic Constitutive Model/Stiffness Property Formulation 

Figure 4.1 shows a schematic of a representative volume element of the SMAHC 
laminate, within the plane of the structure. The principal material coordinates for the 
composite matrix are 1, and 2. Note that in general the figure can represent embedding 
between or within layers and no limitation is placed on SMA orientation relative to the 
composite matrix fiber directions. The following constitutive relation is the fundamental 
equation of the ECTEM developed by Turner “ for the SMA element in the 1 -direction. 


a 


u 



(4.1) 


where, E a is the Young’s modulus of the SMA, S] is the strain in the 1-direction, ai a is the 
coefficient of thermal expansion (CTE). This equation represents the stress-strain relation 
specifically for the SMA material only (hence the subscript a for actuator). The key feature of 

T 

this constitutive relation is the - E a (T )J a Uj ( r]dr term. This tenn embodies the effects of both 

To 

the thermal strain and transformational strain of the SMA. In the constrained recovery 
application, this term is related to the recovery stress, and elastic modulus of the SMA i.e. 
a r (T,s p ), and E a (T), still capturing the nonlinear effects when the temperature is above 
austenitic start ( T>A S ). 


-E a {T)\a u (r)dT = (j r 

To 


or 


T 


To 


°,{T,s p ) 

E a {T ) 


(4.2) 


When the temperature is below austenitic start ( T<A , v ) the thennoelastic relation remains the 
same, and the ECTE of the SMA is due to thermal expansion only, and can be measured 
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experimentally. Therefore, the temperature dependent constitutive relation in the 1 -direction 
for the SMA actuator specific to constrained recovery can be expressed as: 


°u=E a (T) 


*1 - 


i 

\a Xa {v)<E 


T<A s 


To 


(4.3) 
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la 


= E a {T) 


£ x + 


E a (T) _ 


T >A S 


Only the engineering properties of ai a (x) for T < A s and a r , E a for T > A s need to be measured 
in experimental set-ups that imitate the application (pre-strain s p =4%, and boundary 
conditions=clamped) to implement the model. A similar constitutive relation evolves for 
principle material direction-2 (transverse). 


cr 


2 a 



(4.4) 


In this case the transverse CTE 0 C 2 a (x) is not related to the recovery stress, a r , and elastic 
modulus Ea(T), but is still nonlinear due to the changing of the martensite and austenite phases. 

Making the same assumptions, the same stress-strain relationship applies to the 
composite matrix material. Note, the 1,2-CTE’s of the composite matrix are temperature 
dependent. 


T 


*1 " f«l m( T ) dT 

and a 2m = E 2m 

To 




(4.5, 4.6) 


Therefore, the elastic modulus, and CTE in both directions must be measured experimentally 
as functions of temperature. 

Now, the classical mechanics of materials approach to stiffness 38 can be implemented 
to define the material properties of the whole SMAHC lamina in terms of its constituents’ 
measured properties. Imagine that there is only stress present in the 1 direction. The 
assumption is made that the strains in the SMA fiber direction are the same, as the fibers in the 
composite matrix, thus making the plane lamina sections perpendicular to the 1 -direction 
remain plane after stressing. The resultant force in the 1 -direction can be defined as 

°l A l=°la A a+°lm A m ( 4 - 7 ) 
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where, Ai, A a , and A m are the cross-sectional areas of the SMAHC lamina element, SMA 
element, and the composite matrix element, respectively. Then the definition of the average 
stress for the SMAHC lamina element is 

= °la V a + ( 4 ' 8 ) 


where, o a =A a /A i and u m =A m /Ai are the volume fractions of SMA and composite matrix, 
respectively. Then, the stress-strain relationship for each constituent is substituted resulting in: 
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s x - J a x (r)dz 


= En„ 


£ x -\a la {v)di 


+ E x m v n 


*i ~ ja lm (r)dj 


(4.9) 
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To 


To 


From this expression the definition of the SMAHC lamina’s elastic modulus and CTE in the 1- 
direction are formulated. 


i 

E X = E cPa + E mV m and \a x {z)dz 


i i 

E a V a j «1 a + E l m °m J m i T ) dT 


To 


To 


E a U a + E Un V m 


(4.10,4.11) 


where, sgn(« lu ) = 


+ 1 
-1 


T<A S 

T>A S 


This fonnulation for the elastic modulus Ej is physically analogous to the stiffness 
calculation of two springs in parallel. The similar analogy can be made for the modulus 
definition A? in the transverse direction of the SMAHC lamina. In this case, assume a uniform 
stress in the 2-direction only. The strain differs in each constituent, therefore an expression 
representing the total elongation is written. 

ei=e 2m v m +e 2a v a (4.12) 


Substituting for the strain, from equations (2.4, 2.6) results in: 


cr. 


+ J a 2 (r)d t = + j [a 2a (r)u a + a 2m (r)u m r 


J 2 To 


(4.13) 


2m To 


Thus, the Young’s modulus and CTE for the SMAHC lamina in the 2-direction result: 

E E ^ ^ 

E 2 = 77 , ^7 t and J« 2 (r>/r = \[a 2a (r)u a +a 2m {r)u m }tT (4.14) 

\~2m^a E aPm ) To To 
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In this case, the formulation for the Young’s Modulus in the 2-direction (Ei) is analogous to 
two springs in series. The same mechanics of materials approach can be used to define the 
SMAHC lamina’s Poisson’s ratio, V 12 , and shear modulus of a lamina, G 12 . 

The Poisson ratio is similar to the derivation of the Young’s modulus Ej. Assuming 
that there is a unifonn stress in the 1 -direction only on the SMAHC lamina, there is a negative 
defonnation present in the 2-direction, Aw, due to the strain in the 1 -direction. Therefore this 
defonnation in the 2-direction can be represented as 

A w =-Ws 2 (4.15) 


where W is the width of the SMAHC lamina (Fig 4.1). The Poisson’s ratio is defined as 
£ 

v l2 = — - , therefore the transverse defonnation can be represented as 
e \ 

A w =Wv 12 e 1 (4.16) 


Additionally the same transverse defonnation can be represented by the sum of the transverse 
defonnations in the SMAHC lamina constituents, A mW and A aW . 

Wv l2 *l= A mW + A a W ( 4 - 17 ) 


Similar in the analysis for the transverse Young’s modulus E 2 , the deformations A mW and A aW 
are approximately 

^ m w= Wv m v m £\ (4-18) 


A aW = WO a V a £ l 

Thus, upon substitution of (2.18) into (2.17) and division by Wsi, the expression for the 
SMAHC lamina Poisson’s ratio is 

Vn=»mV m +VaVa ( 4 - 19 ) 

Finally, the formulation of the in-plane shear modulus, Gn, of the SMAHC lamina is 
formulated by the assumption that the shearing stresses on the SMA fiber and the matrix are 
the same. The basic assumption can be made that, 

r m = -7T ~ and Ya= -pr ( 4 - 2 °) 

Y2 m Y a 

Let A represent the total shearing defonnation due to shear stress, then 


A = yW (4.21) 

And similarly to the previous fonnulations, the shearing deformation can also be represented 
by the SMAHC lamina constituent’s defonnations. 

A = A m +A a =o m Wy m +o a Wy a (4.22) 
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Therefore, 


y = v y +o y 

/ m/ m a/ a 


(4.23) 


Realization that y = , and substitution of equation (4.20) into (4.23) yields the expression 

G x 2 

for the in-plane shear modulus of the SMAHC lamina. 


G a = 


G„ G 

12m a 

v m G a +o a G nm 


(4.24) 


Now the engineering properties of the SMAHC lamina have all been formulated in 
terms of the individual material constituents (Glass-Epoxy Matrix, and NiTi-SMA), and 
include the thenno-elastic effects due to the ECTE. These properties of the constitutive 
materials can all be measured. Implementing the stress-strain relationships of an orthotropic 
material lamina under plane stress will complete the development of the constitutive relations 
of a SMAHC lamina. 


4.3 Stress-Strain Relation for Orthotropic Material 

In a lamina if there are two orthogonal planes of material property symmetry for a 
material, and symmetry exist relative to a third mutually orthogonal plane, the stress-strain 
relations in coordinates aligned with principal material directions are defined as an orthotropic 
material. If an orthotropic material is subjected to a plane stress state (c3=0, T23=0, X3i=0), the 
three-dimensional strain-stress relation below 


Ai 

D\2 

D n 

0 

0 

0 


^22 

^23 

0 

0 

0 

^13 

^23 

O 33 

0 

0 

0 

0 

0 

0 

D 44 

0 

0 

0 

0 

0 

0 

As 

0 

0 

0 

0 

0 

0 

A, 


(4.25) 


where [d] is the compliance matrix, is reduced to the following stress-strain relation upon 
inversion. 
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where \()\ is the reduced stiffness matrix of a composite lamina. This comes from Hooke’s 
Law in its most basic fonn. This is then implemented for a SMAHC lamina, which includes 
the constitutive relations fonnulated previously in section 4.2. 
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The elements of the reduced stiffness matrix \q\ are related to the engineering properties 
through the following equations. 
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Recall, that these engineering properties of the SMAHC lamina have been defined in terms of 
the lamina material constituents (glass-epoxy, SMA). 
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The above SMAHC constitutive relation (4.27) is for a lamina. This lamina will now 
be considered a general k-th orthotropic layer of a laminate, with an arbitrary orientation angle 
9, from the principal material direction- 1, from Figure 4.1. Below are the transformation 
equations for expressing stress and strain in an x-y coordinate system in terms of stress and 
strain in a 1-2 principal material direction. 
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where [T ] = 


cos 2 # sin 2 # 2sin#cos# 
sin 2 # cos 2 # -2sin#cos# 
-sin#cos# sin#cos# cos 2 # -sin 2 # 


Using the above transfonnations 


(2.30), the thennoelastic constitutive equation for the k th layer of a SMAHC laminate in an 
arbitrary orientation becomes 
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There are some attractive characteristics of this model, which will be mentioned. In 
equation 4.27 it can be seen that the effective CTE in the 1 -direction has an effect on the stress 
in the transverse direction, and vice-versa for the SMA actuator. This indicates that above A s 
temperature the effects of the recovery stress in the 2-direction are included. Also, this model 
not only includes the nonlinear effects of the CTE in the 1 -direction, but also in the 2-direction 
for the SMA actuator. Thus, the behavior of the SMA actuator element in the composite is 
accurately represented in the ECTEM. 

Equation 4.3 1 represents the constitutive relation of the ECTEM in its final form for a 
SMAHC lamina layer. It is evident from the fonnulation that making use of the appropriate 
temperature dependent engineering property measurements and defining a composite laminate 
out of those material constituents is all that is needed for implementation of this model. Any 
commercial finite element code that is capable of analyzing composite laminates with 
temperature dependent material properties can then make use of the ECTEM. Work has been 
done in which an ECTEM research code was compared to ABAQUS, and MSC.Nastran and 
showed excellent correlation 41 . 

In this work, the commercial finite element code, ABAQUS, is used in this way to 
model experimental results of SMAHC beam specimens with dimensions of 
0.5588x0.0254x0.0019 m (22x1x0.078 in). A finite element mesh of 36x4 elements is defined 
with the two center element widths representing the portion of the laminate containing SMA 
material, and the two outer edge element widths representing the portion of the laminate 
containing glass-epoxy matrix material only. Each element is assigned the appropriate material 
stacking sequence and corresponding temperature dependent material properties. This is all 
accomplished in a single input file, in which a number of dynamic/static response analyses can 
be perfonned to accurately represent the experimental conditions. The corresponding material 
characterization, for the material constituents input into the modeling, was previously 
perfonned at NASA LaRC . These material property measurements used for the numerical 
modeling in ABAQUS can be seen in Tables 4.1, 4.2, and 4.3. The Nitinol characterization 
stress-strain curves for 0% and 4% pre-strain, recovery force measurements as a function of 
thermal cycle at 4% pre-strain, and resulting elastic modulus as a function of temperature are 
shown in Figures 4.2, 4.3, 4.4, and 4.5 respectively. These ABAQUS input files will be more 
thoroughly described in chapter 5. 
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Chapter 5 


NUMERICAL MODELING OF SMAHC BEAM SPECIMENS IN THE 
COMMERCIAL FINITE ELEMENT CODE ABAQUS 


5.1 Background 


As previously mentioned, ABAQUS is a commercial finite element code capable of 
analysis of laminated composite structures with temperature dependent material properties. 
Recall that the beam specimen lamination involves the Nitinol ribbon material replacing the 
glass-epoxy matrix over approximately 0.45 inches of the beam width in the 0° layers. This 
can be conveniently modeled by using shell elements of two types; one having glass-epoxy 
layers only, and the other having alternating layers of 100% glass-epoxy and 100% Nitinol 
material. Therefore the ECTE model can easily be implemented by formulating the effective 
composite properties of each constituent material. This type of analysis has a variety of 
element types and solution techniques available to provide sufficient freedom to accurately 
model a wide range of experimental configurations. In this project an initial commercial input 
file for ABAQUS is generated which initially defines the SMAHC structures corresponding 
element, laminate, and material property definitions. Therefore, only a single temperature 
dependent material property data set for each constituent is defined, and used in the correct 
stacking sequence to define the laminate. Thus the laminated SMAHC materials, which are 
composed of orthotropic layers, are defined. The resulting input files in ABAQUS define the 
laminate by the * SHELL SECTION option, and the temperature dependent material properties 
are defined by the *MATERIAL, *ELASTIC, ^EXPANSION options for each material 
constituent (SMA, glass-epoxy matrix). The user is referred to the ABAQUS User’s and 
Reference Manuals and previous work done at NASA LaRC 41 . Once the element and material 
definitions are defined, the ABAQUS input files specify thermal and dynamic loads, and 
solution techniques. The ABAQUS input files will be described briefly below, and a specific 
input file representative of the presented modeling cases can be seen in the Appendix D. 

The unsupported dimensions of the SMAHC beam are 0.4572x0.0254x0.0019 m 
(18x1x0.078 in). This section of the SMAHC beam is modeled by 36x4 type S4 quadrilateral 
shell elements in ABAQUS. As seen in Figure 5.1, the two center element widths (black) 
represent the hybrid composite stacking sequence (45/0/-45/90)2 S with the SMA comprising the 
0° layers, and the two edge element widths (gold) represent the composite stacking sequence 
( 45 / 0 /- 45 / 90 ) 2 s consisting of all glass-epoxy material. Each glass-epoxy layer has a nominal 
thickness of 1.27e-4 m (.005 inches), and each Nitinol ribbon section has a thickness of 1.5e-4 
m (.006 inches). The graphical representation of the beam shows the nodes at each end under 
clamped conditions. Each defined node is represented by the red circle symbol. 

The temperature dependent material properties for each constituent can be seen in the 
ABAQUS input file as seen in the Appendix D. These material properties are taken from the 
material characterization measurements listed in Tables 4.1, and 4.2. It should be mentioned 
that ABAQUS uses secant CTE values as input for thennal expansion and not thennal strain 
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cycle number for the SMA material above the A s temperature. Therefore the appropriate values 
of recovery stress and elastic modulus measurements in Table 4.2 were used to calculate the 
appropriate secant CTE for temperatures above A s temperature. Another important detail worth 
noting is that these material characterization measurements were processed previously at 
NASA LaRC for low thermal cycles 1-4 only. This will be discussed in the subsequent 
sections and impact the comparison of experimental results to the numerical simulations. The 
appropriate CTE values for the Nitinol ribbons below the A s temperature were taken from 
published data as were the transverse direction CTE values for martensitic temperature 6.6e-6 
/°C (3.67e-6 /°F) and austenitic temperature lle-6 /°C (6.11e-6 /°F) . The transverse CTE 
values have been linearly interpolated between these two points to represent the values through 
the temperature range where the SME occurs. 

Once the geometry, finite element mesh, and material properties are defined, the 
physical conditions of each experiment can be simulated. Some of the experimental 
measurements will serve as input into the numerical simulations. The ABAQUS input files 
consist of a sequence of several analysis steps progressing through the application of the 
thermal load. At certain load levels a desired equilibrium state is achieved, where the random 
response analysis is simulated for those states, thus determining the dynamic response. In each 
individual input file the actual measured thennal load from the thennal imaging camera for 
each thermal step to 121.1 °C (250 °F) is applied to the SMAHC beam with clamped boundary 
conditions. Throughout the application of the thermal load a dynamic inertial loading is 
defined that represents the .25g RMS base excitation used in the dynamic experiments. The 
actual base excitation PSD measured from each experiment, see Figure 3.3 for an example, is 
input into the analysis by the *PSD-DEFINITION in ABAQUS. During the linear random 
response step analysis at the specified equilibrium states a modal approach is used 
incorporating the first 10 modes, in order to account for residual effects. The modal parameter 
estimates for mode 1 and 3 measured during the dynamic response experiments are used as 
input for damping estimates through modes 1 through 10 using the *MODAL DAMPING 
option. Additionally, the initial geometric imperfection measured by PMI measurements for 
the full-field displacement of the beam is input into the analysis by the ^IMPERFECTION 
option, to assist in the nonlinear static solution procedure of the buckling/post-buckling of the 
beam. Recall the PMI measurements were for the centerline of the beam length only, but are 
representative of the full-width of the beam. Therefore a MATLAB m-file was written which 
linearly interpolates the beam length centerline initial imperfection measurements to nearest 
nodes corresponding to the full field of the SMAHC beam (See m-file and corresponding data 
set in Appendix C). 

In the subsequent sections the numerical simulation input files will contain the 
previously discussed experimental measurements, which are beam specimen and run specific. 
This should improve the accuracy of the numerical simulations. Results for the analysis of 
nonlinear static, and linear random dynamic responses will be presented and compared to 
experimental results. 

5.2 Static Response Modeling Results in ABAQUS 

The plots of midspan deflection of the SMAHC beam vs. average centerline 
temperature are presented for the experimental PMI measurements and the numerical 
predictions in ABAQUS. Recall, from the static thermoelastic response experiments that the 
thermal buckling maximum midspan deflection varies greatly from run to run, and may be 
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attributed to moisture absorption, and varying pre-stress. Thus the specifically chosen 
numerical simulations will be given the initial imperfection data for runs that correspond to the 
middle of the varying deflection range where repeated experimental results occurred. 

As previously mentioned, during the numerical simulations the specific material 
property measurement data for each constituent (SMA, Glass-Epoxy) is input into the model. 
Looking at figure 4.4 the recovery force data shows that A s temperature begins when stress is 
starting to be generated at 28 °C (82.4 °F). Isolength (constant strain) tests were performed on 
a nominal .254 m (10 inch) length of Nitinol ribbon material to obtain these recovery stress 
results as a function of pre-strain, and thermal cycle. Thermocouple’s were used to measure 
the corresponding temperature. There is a large error associated with thennocouple 
measurement at high temperatures due to ‘fin’ effect. There is also error due to ‘fin’ effect 
even at low temperature if the thennocouple is attached to a small surface area, as the case with 
the SMA ribbon. This recovery stress data for thennal cycle 4 is used for numerical modeling, 
but it is suspected that there is some error in the conesponding temperatures. This is validated 
by a discrepancy in the measured A s temperature by differential scanning calorimetry (DSC) 
measurements. 

Figure 4.6 shows the DSC results for small SMA ribbon samples (<10mg). The 
exothermic and endothermic peaks conespond to the transfonnation of the SMA material. 
During the first-heat there is an ill-behaved endothermic peak associated with the austenitic 
transformation recovery shear defonnation induced from cutting the samples. This peak was 
not repeatable after the first heat, where as the other endothennic peak was. Therefore the 
measured A s temperature is near 40 °C (105 °F). These measurements are accurate and indicate 
a possible error with the recovery stress measurements. This will be shown in the results of the 
modeling. 

For each static thennoelastic response simulation the corresponding experimental 
measurement for a given initial imperfection, and thennal load is plotted. The numerical 
simulations for the original material property measurements, and after shifting the recovery 
force data 10-20 °F to account for the possible measurement enor are included. Material 
property data is input into ABAQUS as elastic modulus data (E) and coefficient of thermal 
expansion (CTE) data. Recall, from the model, that the Nitinol effective CTE values are 
obtained from the measured E and the measured recovery stress (a r ) at temperatures above A s . 
Therefore the a r data is shifted and the conesponding CTE values are obtained from equation 
4.2. These shifted CTE values are then used as input into ABAQUS to adjust for the 
measurement enor. 

In Figure 5.2 the predicted static response for Beam 2 Run 6 is shown. This numerical 
prediction for clamped boundary conditions, and the original recovery stress (CTE) 
measurements indicates that the SMA recovery force is overcoming the compressive forces 
due to expansion at a lower temperature than indicated by experiment, thus resulting in a lower 
max deflection, and smaller post-buckled temperature range (see green dashed line). The 
recovery stress measurements were shifted 1 5 °F and the corresponding new CTE values were 
input instead. This results show great qualitative agreement between the numerical simulation 
and the experimental measurement (See blue dashed line). The same results were produced in 
multiple simulations for Beam 2 and Beam 3 (Figures 5.3 - 5.7). Some cases required varying 
shifts between 10-20 °F. These results indicate that the pre-stress initially present could be 
varying for run to run slightly affecting the recovery force behavior. These static experimental 
results are based on thennal cycles 5-14 for beam 3 and cycles 1-9 on beam 2. The recovery 
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force measurement data was only captured for thermal cycles 1-4. Therefore thennal cycle 4 
data was used in modeling. This could be attributed to the small differences between 
experiment, and simulation. The recovery force data adjustment seems to give great qualitative 
agreement with experiment and indicates the sensitivity of the buckling behavior to the 
interaction between the SMA recovery stress generated and the compressive forces due to 
thermal expansion. 

The recovery stress temperature shift now activates the recovery forces at the correct 
corresponding temperature to more accurately capture the actuation, which returns the midspan 
deflection to zero. In all simulations the measurements of the initial imperfection, the thermal 
load, and material property measurements were input directly into the model to achieve 
accuracy. These results verify the capability of the commercial finite element code ABAQUS, 
to accurately predict the static thennoelastic response of SMAHC structures by implementing 
the ECTE model. Additional material characterization measurements could be made for the 
recovery stress of the constrained Nitinol actuator ribbons using infrared thennal imaging to 
produce even more accurate material property data to be used in the model. This could even 
further improve the accuracy of the numerical prediction with the experimental results. 

5.3 Dynamic Response Modeling Results in ABAQUS 

Experimentally measured and numerically predicted displacement PSD’s for 
measurement position 2, representing the midspan accelerometer, are shown for various 
temperature values in Figures 5.8-5.15. Results are presented only for beam 2 run 8 and beam 
3 run 8, but are representative of all other dynamic results for that measurement position, and 
beam specimen. Dynamic response PSD’s are presented for both beam specimens at 21.7 °C 
(71 °F), 71.1 °C (160 °F), 93.3 °C (200 °F), and 121.1 °C (250 °F). These temperature values 
demonstrate the stiffening of the beam due to the SMA actuation and the ability of the 
numerical simulation to accurately capture the effect. For these numerical simulations the 
adjusted recovery stress material data was input for the material properties (15 °F shift). 

Recall, from experimental results in chapter 2 it was discovered that the boundary 
clamped G-10 insulation was unsymmetrical and the boundary conditions are contributing to 
the initial imperfection of the beam. Also the experimental initial imperfection measurements 
show that the boundary conditions are dominating in the initial shape of the beam. The fact 
that the G-10 insulation thickness was not unifonn for all samples indicates that the boundary 
conditions are most likely not exhibiting perfectly clamped conditions, and would be difficult 
to accurately model. The unsymmetrical G-10 insulation clamps likely allow some torsional 
rotation of the beam at the boundary. This was attempted to be accounted for in the model by 
applying torsional spring elements to the 5th degree of freedom (out of plane of the beam) of 
the end nodes at the boundaries of the beam model. This is defined in ABAQUS by the 
^SPRING function. Torsional springs with a rate of 60 (in-lbf/rad) were applied to match the 
fundamental frequency. 

As indicated in previous work 18 , in the physical system the beam continues through the 
mechanical grips and the Nitinol leads are terminated beyond that in the electrical grips. 
Therefore the main restraint for the Nitinol is the mechanical tennination at the electrical grips. 
Thus, the compliance of the entire support structure needs to be considered. Although, Garolite 
G-10 material was used as insulation between the electrical and mechanical grips, it is likely 
that there is still some compression occurring, partially relieving the recovery stress. This effect 
was captured in the model by applying axial springs with a spring rate of 7000 (lbf/in) to each 
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of the nodes at the boundary of the beam model in degree of freedom 1 (in-plane). This spring 
rate was obtained by matching the fundamental frequency of the simulation to the measured 
fundamental frequency at the elevated temperature of 250 °F. 

The same explanations described in the previous static thennoelastic response modeling 
sections correspond to the over prediction of stiffness initially seen in dynamic results. The 
application of the boundary springs previously discussed produce the PSD comparison in 
Figures 5.8-5.15, for beam 2 run 8 and beam 3 run 8 at position 2 (midspan). The displacement 
PSD predictions show excellent agreement with the experimental results. It can be seen that at 
21.7 °C (71 °F), 71.1 °C (160 °F), 93.3 °C (200 °F), and 121.1 °C (250 °F) temperatures, the 
ABAQUS predicted responses closely correspond with the experimental modal frequency 
values and RMS displacement values. In Figures 5.16 and 5.17, the same temperature values 
are compared to the experimental measured PSD’s for positions 1 and 3, along the length of the 
beam all in one plot. Again the ABAQUS predictions show excellent agreement with the 
experimental results in modal frequency values, and RMS displacement values. 

In Figures 5.18-5.23 displacement PSD’s of the numerical simulations and 
experimental dynamic responses are compared for temperatures corresponding to the buckled 
and post-buckled range of the SMAHC beam specimen. These particular data sets correspond 
to beam 3 run 1 of the dynamic experiments. This is the only dynamic experimental run in 
which dynamic data was captured through this temperature range of buckling/post-buckling. 
The static deflection behavior of the beam doesn’t correspond closely the numerical simulation 
of the static deflection for this particular run in that temperature range. Therefore the PSD 
comparison at these temperatures is expected to suffer as well. If dynamic data was captured 
for this temperature range for the thermal cycles presented in Figures 5. 2-5. 7 then the post- 
buckled temperature range PSD’s would be expected to compare reasonably well. Through 
this temperature range the beam goes through the buckling phenomenon in which the stiffness 
theoretically drops to zero. Once the beam is buckled the stiffness is increased, until the SMA 
recovery forces begin to return the buckling deflection back to a flat configuration. At this 
point the stiffness of the beam is decreasing. So during this temperature range the SMAHC 
beam exhibits complex dynamic behavior, and would be difficult to model accurately. The 
PSD comparison in Figures 5.18-5.23 demonstrates the difficulty in modeling this behavior 
with the current analysis, and the complex dynamic behavior evident in the experimental 
PSD’s can be seen. 

In Figure 5.18 the displacement PSD at ambient temperature 21.7 °C (71 °F) is plotted 
for position 3 along the length of beam 3. This particular dynamic response demonstrates some 
nonsymmetrical modes as seen by the additional spectral peaks in the experimental PSD. Also 
nonlinearity is demonstrated by peak shifting and peak broadening in the experimental PSD’s. 
ABAQUS numerical model shows good agreement with the stiffness prediction even though a 
linear dynamic response is predicted. Recall that a linear random analysis is being 
implemented in ABAQUS to predict the dynamic response. That explains the difference in the 
PSD’s comparison, but the modal frequency peaks still match indicating a good prediction of 
the stiffness state of the SMAHC beam. In Figure 5.19 the displacement PSD comparison for a 
temperature of 37.7 °C (100 °F) is shown for Beam 3 Runl at position 3. At this temperature 
the beam has buckled. The experimental PSD’s show again the complex nonlinear behavior of 
the dynamic response, and shows a slight increase in the fundamental frequency from the 
ambient state. The numerical simulation in ABAQUS does not agree with the experimental 
results in that it is predicting a large decrease in stiffness by the decrease in the modal 
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frequency values. The same prediction error is shown in Figure 5.20 at a temperature of 43.3 
°C (130 °F). At this temperature, recovery stresses are being generated because it is above A s 
temperature, but they have not overcome the compressive forces due to thermal expansion, 
thus the beam is still in a buckled state. The experimental PSD is again showing complex 
dynamic behavior by the nonlinear peaks. In Figure 5.21 the PSD comparison at 60 °C (140 
°F) shows a good agreement in modal frequency values. At this temperature state the beam has 
been just been rendered flat by the SMA recovery forces. Since the beam has just been 
returned flat it is still demonstrating some nonlinear dynamic behavior, but ABAQUS is able to 
make better predictions now using its linear analysis. Then in Figure 5.22 at 65.55 °C (150 °C) 
the experimental PSD shows that the linear dynamic behavior has returned as the beam is flat 
and an increase in stiffness is occurring. The ABAQUS numerical simulation now is in good 
agreement with the behavior for the rest of the temperature range and accurately predicts the 
stiffness of the structure as seen in subsequent Figure 5.23 at 71.1 °C (160 °F). 

So the dynamic response modeling results demonstrate the ability of ABAQUS to 
accurately model the dynamic behavior of the SMAHC beam specimens’ experimental results 
by using a linear random analysis. The adaptive stiffening effect and RMS values of the 
vibration are accurately predicted by input of the measured modal parameters, the random 
excitation from experiment, and the measured material properties of the constituents directly 
into the numerical model. Accurate simulations are produced in the temperature range below 
and above the buckling range, and demonstrate the stiffening effect of the embedded SMA 
actuators. An attempt to model the dynamic behavior through the temperature range, where 
buckling occurs, demonstrates that there is complex dynamic behavior occurring when the 
SMAHC beam is in this state. If the static deflection behavior of the beam corresponds closely 
to the numerical simulations as in Figures 5.2-5. 7, it expected that the numerical PSD 
prediction would correspond closely in this post-buckled temperature range except at the start 
of the buckling and at the point of returning flat. These two locations are theoretically the only 
cases where linear analysis fails due to the nonlinear and complicated dynamics. This limits 
the accuracy of a linear random analysis to predict the behavior. It would be desirable perhaps 
for application advantages to use Nitinol actuator ribbons with a Ni and Ti composition which 
produces an As temperature below ambient temperature to avoid the buckling behavior of the 
structure. 
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Chapter 6 


CONCLUSIONS 


This study consisted of a thorough experimental investigation into the response 
behavior of SMAHC beams. The attractive functionality of SMA’s was demonstrated as an 
adaptive stiffening mechanism to composite beam specimens in all experimental cases. The 
ECTE constitutive model was implemented into the commercial finite element code ABAQUS, 
to accurately predict the response of the SMAHC structure. The experimental measurement 
investigations uncovered important contributing factors to the behavior of the SMAHC beam, 
such as initial geometric imperfections, contributions of the boundary conditions, distributed 
thermal load over the beam, and the real base excitation that was used to induce the dynamic 
response. These measurements, and discoveries were able to be accounted for in the numerical 
modeling to produce an accurate prediction of the SMAHC beam responses, and further 
validate the use of the ECTE model to be used in the commercial finite element code 
ABAQUS. 

Numerous experimental tests were perfonned to detennine the thermoelastic response 
in a clamped boundary configuration with a thermal load applied via resistive heating. The out 
of plane deflection of the midspan was measured as a function of temperature and uncovered a 
possible dependency of maximum deflection on moisture absorption in the hybrid composite 
structure. Excellent correlation was demonstrated between the PMI measurement technique 
and the laser displacement transducer measurements. The PMI measurements enabled the 
discovery of the initial geometric imperfections of the beam, and the extreme involvement of 
the beam fixture boundary conditions on this imperfection. This validates the accuracy of the 
PMI measurement technique for this application, and gives the valuable understanding of the 
full-field displacement of the structure during testing. 

The SMAHC beam specimen’s dynamic response was investigated with a generated 
thermal load by resistive heating, and application of a 10-400 Hz random base excitation. The 
dynamic measurements were post-processed into FRF’s, PSD’s, RMS values, and coherence 
functions. The post-processing of the experimental results not only demonstrated the adaptive 
stiffening mechanism of the SMAHC beams, but also enabled the modal parameter estimation 
of damping, and modal frequency of each mode. The displacement RMS values were shown 
to decrease with temperature. The modal damping estimates decreased once the Nitinol ribbon 
was above A s temperature, until higher temperatures near 121.1 °C (250 °F) were reached, 
where a slight increase in damping estimates occurred most likely due to the softening of the 
heated glass epoxy material. 

The precise experimental measurements were used as input into the commercial finite 
element code ABAQUS to accurately predict the static and dynamic responses of the SMAHC 
beam specimens. In the static thennoelastic response comparisons of numerical prediction and 
experimental results, ABAQUS showed excellent qualitative agreement with experiment for 
multiple runs. It was discovered that there is most likely a temperature correlation error in the 
Nitinol recovery force measurements of 10 - 20 °F. When this error is accounted for in the 
CTE material property input values excellent agreement with the experimental results is 
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demonstrated. This indicates that it is critical that the appropriate material characterization 
measurements are accurately measured for the state that the embedded SMA material is in. 

The numerical prediction of the dynamic response of the SMAHC beams proved to be 
in excellent agreement with experimental measurements after the application of boundary 
springs, and shifting of the CTE material property values. The displacement PSD’s 
comparison showed agreement for all positions 1, 2, and 3 along the length of the beam for 
both beam specimens. The results were achieved by a linear random analysis in ABAQUS 
using accurately measured modal parameters, material property data, thennal loads, and 
random base excitation as input into the model. An additional investigation into the dynamic 
behavior of the beam at temperatures corresponding to the buckling phenomena (80 °F - 140 
°F) showed the limitations of the numerical model using a linear random analysis to capture the 
complex dynamic behavior which occurs through this temperature range. 

The motivation for the static and dynamic response investigation of SMAHC beams is 
to demonstrate the decrease in the dynamic response amplitude when under random 
excitations. Random excitation of a structure can lead to sonic fatigue and undesirable noise 
transmission characteristics of a structure. By demonstrating a reduction in the response, and 
the adaptive stiffening mechanism introduced by the SMA actuators, the implications of 
vibration and noise transmission reduction are evident. It is possible that finding an optimal 
lay-up or embedding position of the SMA actuators to demonstrate most efficient control over 
the structure could result in increased performance of the SMAHC structure. Also, the 
fabrication of other SMAHC structures such as a panel specimen was previously done at 
NASA LaRC, although experimental validation wasn’t completed. Numerical studies have 
shown that a weight-efficient SMAHC panel demonstrates superior dynamic response and 
sonic fatigue behavior over other conventionally noise-treated panels. Also, the SMAHC panel 
numerically shows superior transmission loss (TL) characteristics in the low frequency, 
stiffness-controlled region. The weight-efficient stiffening effect has significant implications 
for dynamic response and sonic fatigue reduction. For further research, the fabrication of 
another SMAHC panel could be done, in which an extensive experimental investigation could 
be perfonned measuring the vibration response, acoustic radiation, and transmission loss 
characteristics using the more thorough measuring techniques used in this study. The acoustic 
radiation and transmission loss experimental investigation could produce a multitude of 
additional discoveries. With the addition of increased commercial code capabilities the 
numerical simulation study could also be investigated to compare. 

In conclusion, this study uses fundamental engineering property measurement to 
accurately model SMAHC beam specimens. The experimental investigation of the static and 
dynamic response of the SMAHC structures was performed using accurate measuring 
techniques, which uncovered interesting and complex behavior of the SMAHC beam 
specimens. The SMA actuators demonstrated an adaptive stiffening effect to the structure, 
which was successfully, modeled in the commercial finite element code ABAQUS. 
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Chapter 8 


TABLES 


Table 3.1- Beam 2 Run 1 dynamic response RMS displacement values as a function of 
temperature for positions 1, 2, and 3 along the length of the beam. 



Post 

Pos2 

Pos3 

Temp 

RMS 

RMS 

RMS 

°F 

Displ., in 

Disp, in 

Displ., in 

73 

.838e-2 

1.30e-2 

,645e-2 

150 

.656e-2 

.985e-2 

.586e-2 

160 

.523e-2 

.757e-2 

.477e-2 

170 

,441e-2 

.627e-2 

.404e-2 

180 

.473e-2 

.607e-2 

.435e-2 

190 

.343e-2 

,482e-2 

.317e-2 

200 

.320e-2 

.701e-2 

.297e-2 

210 

.297e-2 

.415e-2 

.275e-2 

220 

.267e-2 

.373e-2 

.247e-2 

230 

.260e-2 

.364e-2 

,241e-2 

240 

.254e-2 

.355e-2 

.235e-2 

250 

.254e-2 

.356e-2 

.235e-2 


38 




Table 3.2- Beam 3 Run 2 dynamic response RMS displacement values as a function of 
temperature for positions 1, 2, and 3 along the length of the beam. 



Post 

Pos2 

Pos3 

Temp 

RMS 

RMS 

RMS 

°F 

Displ., in 

Disp, in 

Displ., in 

73 

,941e-2 

1.69e-2 

.812e-2 

150 

.652e-2 

1.03e-2 

.669e-2 

160 

.376e-2 

.573e-2 

.336e-2 

170 

.366e-2 

.551e-2 

.389e-2 

180 

.394e-2 

.587e-2 

,427e-2 

190 

,343e-2 

.508e-2 

.381e-2 

200 

.331e-2 

.489e-2 

.347e-2 

210 

.294e-2 

,434e-2 

.278e-2 

220 

.263e-2 

.388e-2 

.259e-2 

230 

,251e-2 

.370e-2 

,223e-2 

240 

,243e-2 

.358e-2 

,228e-2 

250 

.239e-2 

.352e-2 

,218e-2 
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Table 3.3- Beam 2 dynamic response modal parameter estimates for runs 1, 2, and 3 



Runl 

Run2 

Run3 

Temp , °F 

fi , Hz 

Ci, % 

f 3 , Hz 

C.3, % 

fi , Hz 

Ci, % 

f 3 ,Hz 

C3, % 

fi,Hz 

Ci, % 

f 3 , Hz 

Cj,% 

73 

30.1 

1.06 

160.2 

.66 

30.9 

.81 

163.3 

.68 

30.5 

.87 

162.1 

.73 

150 

59.6 

.31 

223.2 

.66 

67.1 

.17 

240.0 

.39 

66.1 

.22 

237.8 

.43 

160 

80.0 

.17 

271.9 

.26 

81.7 

.17 

276.1 

.25 

81.9 

.20 

276.7 

.29 

170 

91.3 

.14 

302.0 

.27 

92.3 

.15 

304.5 

.26 

92.9 

.14 

306.4 

.29 

180 

99.5 

.15 

323.7 

.19 

100.8 

.13 

327.3 

.18 

101.5 

.13 

329.1 

.18 

190 

106.9 

.13 

343.9 

.15 

107.5 

.15 

345.5 

.16 

108.1 

.13 

347.1 

.14 

200 

112.9 

.16 

360.2 

.17 

113.1 

.16 

360.7 

.16 

112.6 

.17 

359.6 

.14 

210 

117.2 

.18 

371.9 

.18 

116.9 

.16 

371.2 

.16 

116.1 

.16 

369 

.15 

220 

120.3 

.20 

380.9 

.27 

119.5 

.17 

378.5 

.20 

118.4 

.17 

375.4 

.17 

230 

122.5 

.21 

387.4 

.24 

121.1 

.18 

383.1 

.27 

120.0 

.18 

379.9 

.25 

240 

124.1 

.21 

391.6 

.19 

122.0 

.18 

385.8 

.23 

121.1 

.19 

383.1 

.27 

250 

125.1 

.21 

394.2 

.18 

122.8 

.21 

388 

.22 

121.2 

.19 

384.7 

.25 
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Table 3.4- Beam 3 dynamic response modal parameter estimates for runs 1, 2, and 3 



Runl 

Run2 

Run3 

Temp , °F 

fi , Hz 

Ci, % 

f 3 , Hz 

C.3, % 

fi , Hz 

Ci, % 

f 3 ,Hz 

C3, % 

fi,Hz 

Ci, % 

f 3 , Hz 

Cj,% 

73 

30.3 

.72 

162.6 

.56 

29.9 

.81 

161.9 

.59 

31.5 

.76 

163.9 

.66 

150 

60.5 

.24 

222.4 

.68 

62.8 

.24 

225.2 

.97 

67.3 

.21 

237.9 

.84 

160 

79.8 

.29 

269.4 

.21 

80.7 

.29 

272.3 

.40 

82.7 

.33 

276.5 

.42 

170 

91.3 

.12 

298.9 

.34 

92.5 

.20 

301.9 

.35 

93.5 

.21 

304.1 

.42 

180 

100.5 

.10 

324.5 

.42 

101.3 

.15 

324.6 

.58 

101.9 

.14 

326.7 

.75 

190 

107.1 

.09 

342.5 

.18 

108.2 

.14 

346.4 

1.13 

108.3 

.14 

346.2 

1.18 

200 

112.1 

.11 

356.3 

.14 

113.7 

.16 

362.4 

.38 

113.0 

.15 

360 

.52 

210 

116.5 

.15 

368.5 

.12 

117.8 

.16 

372.2 

.29 

116.8 

.18 

369.8 

.34 

220 

119.1 

.18 

375.6 

.12 

120.3 

.20 

378.7 

.44 

119.4 

.23 

375.9 

.35 

230 

121.1 

.21 

381.2 

.12 

122.2 

.21 

383.8 

.86 

121 

.23 

380 

.74 

240 

122.5 

.24 

384.6 

.17 

123.3 

.22 

388 

1.10 

121.8 

.22 

383.1 

1.12 

250 

123.9 

.25 

388.8 

.28 

123.8 

.20 

389.9 

1.31 

122.3 

.21 

384.9 

.78 
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Table 4.1- Thermomechanical material property measurement of the Fiberite E- 
glass/934 epoxy material system in principal material direction. 


Temperature 

Ei 

E: 

"12 

Chi 

tti 

°2 

°C 

GPa 

GPa 

— 

GPa 


t n/°C 

21.1 

49.30 

20.00 

0.29 

9.65 

5.22 

10.80 

26.7 

49.30 

19 86 

0.29 

9.58 

5.83 

13.81 

37 8 

49.16 

19.44 

0.29 

9.24 

6.95 

21.02 

48.9 

49.02 

18.96 

0.29 

8.89 

7.22 

24.14 

60.0 

48.82 

18.48 

0.29 

8.55 

7.00 

24.66 

65.6 

48.75 

18.20 

0.29 

8.41 

6.80 

24.28 

71.1 

48.75 

17.79 

0.29 

8.27 

6.62 

23.90 

82.2 

48.68 

17.03 

0.29 

7.93 

6.34 

22.79 

93.3 

48.61 

16.20 

0.29 

7.58 

6.25 

22.00 

104.4 

48.61 

15.31 

0.29 

6.76 

6.39 

21.92 

115.6 

48.54 

14.41 

0.29 

6.00 

6.77 

22.66 

121.1 

48.54 

14.00 

0.29 

5.58 

6.97 

23.33 

126.7 

48.61 

13.45 

0.29 

5.17 

7.18 

24.01 

137.S 

48.68 

12.41 

0.29 

4.27 

7.42 

25.47 

148.9 

48.82 

11.38 

0.29 

3.45 

7.63 

26.28 
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Table 4.2- Thermomechanical material property measurements of the Nitinol Ribbon 
actuator material. 



Recover*’ Stress 



Temperature 

Cycle 2 

Cycle 3 

Cycle 4 

E 

02 

°C 

MPa 

MPa 

MPa 

GPa 

,i(P c 

21.1 

0 

0 

0 

27.17 

6.606 

26.7 

2.130 

1.276 

2.130 

24.82 

6.606 

32.2 

7.233 

5.530 

7.233 

22.41 

6.606 

37 S 

16.18 

13.20 

14.05 

20.06 

6.606 

43.3 

31.92 

28.95 

28.52 

25.72 

6.606 

48.9 

70.65 

70.23 

78.74 

31.37 

7.236 

54.4 

135.4 

137.9 

151.9 

36.96 

7.866 

60.0 

192.0 

199.6 

217.9 

42.61 

8.496 

65.6 

240.5 

256.6 

269.4 

48.27 

9.10S 

71.1 

297.1 

301. S 

309 8 

54.S8 

9.738 

76.7 

341.8 

341. S 

341.3 

61.43 

10.37 

82.2 

374.5 

373.7 

372.0 

64.19 

11.00 

87.8 

402.2 

399.2 

399.2 

63.16 

11.00 

93.3 

425.6 

419.2 

4188 

62.06 

11.00 

98.9 

449.9 

438.8 

437.5 

63.92 

11.00 

104.4 

469.0 

451.6 

448.2 

65.78 

11.00 

110.0 

481.0 

464.8 

458.8 

67.64 

11.00 

115.6 

491.6 

475.8 

467.3 

69.50 

1100 

121.1 

500.1 

485.6 

472.9 

71.36 

11.00 

126.7 

508.6 

490.3 

478.8 

70.81 

11.00 

132.2 

515.0 

496.7 

481 8 

70.33 

11.00 

137.8 

520.5 

500.5 

485.2 

69.78 

11.00 

143.3 

523.5 

502.2 

487.3 

69.29 

11.00 

148.9 

527.8 

503.9 

490.7 

6S.74 

11.00 
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Table 4.3- Nitinol Ribbon Transformation Temperatures measured by DSC. 



Temperature 

Transformation Marker 

°C 

°F 

Austenite Start 

45 

113 

Austenite Finish 

60 

140 

Martensite Start 

17 

62.6 

Martensite Finish 

0 

32 
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Chapter 9 


FIGURES 




Martensite (twinned) Martensite (M+) (detwinned) Austenite 

Figure 1.1- SMA Lattice Structure at different phases. 


Heating 



Cooling 


Figure 1.2- SMA body under tensile loading/unloading at low temperature, exhibiting 

the shape memory effect upon heating 4 . 


45 




Temperature 


Figure 1.3- Schematic of martensitic and reverse transformations 5 . 







Figure 1.4- Stress-Strain Curves for Shape Memory Alloys at various intermediate 

temperatures 6 . 
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quasiplastic 


pseudoelastic 


Figure 1.5- Load-Deformation diagrams for shape memory alloy at various 
temperatures (fully martensitic-quasiplastic) and (fully austenitic-pseudoelastic) 4 . 



Figure 1.6- Nitinol recovery stress vs. temperature and initial strain 7 . 



Figure 1.7— Recovery stress vs. temperature, thermal cycle for nitinol ribbon samples at 

4% pre-strain 8 . 



Figure 1.8- Variation of austenite start temperature with weight content (%) of Nickel 

in Nitinol Alloy 10 . 


48 


1 


0.5 


1 " 

r 

I- 22" 

Figure 2.1- 22”xl”x0.080” Glass Epoxy/SMA Hybrid Beam (45/0/-45/90)2 S Lamination 
with SMA embedded in the 0° lamina only 8 . 





Figure 2.2- Consolidated SMAHC structure after curing cycle, and machined SMAHC 

beam specimens used for experiments 8 . 



Figure 2.3- Static Experimental Set-up: SMAHC Beam3 installed in the beam fixture 

(front view) 
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Figure 2.4b- Mechanical and electrical grip of beam fixture (back view). 
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Figure 2.5- Static Thermoelastic Response Experimental Component Schematic. 
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Figure 2.6- Static Thermoelastic Response Experimental Configuration Overview. 


51 












Midspan Displacement (in) 



Figure 2.7- Simple PMI System Schematic 35 . 
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Figure 2.8- SMAHC Beam2 Laser Displacement Experimental Results of Midspan for 

Runs 1-9. 
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Figure 2.9- SMAHC Beam3 Laser Displacement Experimental Results for Midspan 

Runs 1-10. 



Figure 2.10- Experimental Configuration of initial imperfection measurements. 
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Beam Length (in) 


Figure 2.11- Laser Translation versus PMI initial imperfection measurements of 
SMAHC beams installed in beam fixture. 
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Figure 2.12- Laser Translation Measurements of SMAHC Beam3 for three different 
boundary condition configurations of the Garolite G-10 insulation blocks. 
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Figure 2.13- Garolite G-10 Insulation Blocks Configurations for Laser Displacement 
Transducer Experiments (Top view of beam fixture boundaries) 
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Figure 2.15- PMI differential image data and line fit for the centerline out of plane 

deflection for Beam3 Run3 at T=160°F. 
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Figure 2.16- PMI full field imaging of the SMAHC Beam2 Runl 
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Figure 2.17- PMI full field imaging of the SMAHC Beam 2 Run 4 
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Figure 2.18- PMI full field imaging of the SMAHC Beam 2 Run 8 
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Figure 2.19- PMI full field imaging of the SMAHC Beam 3 Run 1 
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Figure 2.20- PMI full field imaging of the SMAHC Beam 3 Run 4 
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Figure 2.21- PMI full field imaging of the SMAHC Beam 3 Run 8 
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Figure 2.22- Beam2 PMI vs. Laser Static Deflection Measurements 
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Figure 2.23- Beam3 PMI vs. Laser Static Deflection Measurements 



Figure 3.1- SMAHC beam dynamic response experimental configuration. 
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Figure 3.2- SMAHC beam specimen installed, with mounted beam response 
transducers and fixture transducers (back view). 



Figure 3.3- Acceleration PSD for the base excitation of the beam with a RMS value of 
.25g (96.5 in/s 2 ). 
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Figure 3.4- SMAHC beam dynamic response experimental component schematic. 
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Figure 3.5- Acceleration vs. Displacement PSD generated by attached m-file, which post- 
processes the acceleration time history data. 



Figure 3.6- SMAHC Beam 2 PSD’s for selected temperatures at position 1. 
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Figure 3.7- SMAHC Beam 2 PSD’s for selected temperatures at position 2 (midspan). 



Figure 3.8- SMAHC Beam 2 PSD’s for selected temperatures at position 3. 
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Figure 3.9- SMAHC Beam 3 PSD’s for selected temperatures at position 1. 
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Figure 3.10- SMAHC Beam 3 PSD’s at selected temperatures for position 2 (midspan). 
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Figure 3.11- SMAHC Beam 3 PSD’s at selected temperatures for position 3. 



Figure 3.12- Modal frequency variation of modes 1, and 3 as a function of temperature. 
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Figure 4.1- Representative volume element of SMAHC lamina 
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Figure 4.2- Applied stress vs. strain for .09x.006xl0 inch sample with 0% pre-strain at 

ambient temperature and 120°C (248 °F) 18 . 
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Figure 4.3- Applied stress vs. strain for .09x.006xl0 inch sample with 4% pre-strain at 

ambient temperature and 120°C (248 °F) 18 . 



Figure 4.4- Recovery Stress vs. Temperature of the Nitinol ribbon (pre-strained 4% 
and released) for thermal cycles 2, 4, and 50 18 . 
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Figure 4.5- Modulus for Nitinol ribbon prestrained 4% and released as a function of 

i o 

temperature . 



Figure 4.6- Differential Scanning Calorimetry (DSC) Results for the Nitinol ribbon 
sample showing transformation temperatures. 
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Figure 5.1- ABAQUS Finite Element Mesh of the SMAHC beam with clamped 

boundary conditions. 



Figure 5.2- Static Response SMAHC Beam 2 Run 6 ABAQUS vs experimental results 
for original recovery force measurements, and 15 °F shift in recovery force input data. 
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Figure 5.3- Static Response SMAHC Beam 2 Run 7 ABAQUS vs. experimental results 
for original recovery force measurements, and 10 °F shift in recovery force input data. 



Figure 5.4- Static Response SMAHC Beam 2 Run 8 ABAQUS vs. experimental results 
for original recovery force measurements, and 15 °F shift in recovery force input data. 
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Figure 5.5- Static Response SMAHC Beam 3 Run 3 ABAQUS vs. experimental results 
for original recovery force measurements, and 15 °F shift in recovery force input data. 



Figure 5.6- Static Response SMAHC Beam 3 Run 6 ABAQUS vs. experimental results 
for original recovery force measurements, and 15 °F shift in recovery force input data. 
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Figure 5.7- Static Response SMAHC Beam 3 Run 8 ABAQUS vs. experimental results 
for original recovery force measurements, and 20 °F shift in recovery force input data 



Figure 5.8- Beam 2 Run 8, position 2 (midspan) PSD comparison at ambient (Af=.25 

Hz). 
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Figure 5.9- Beam 2 run 8, position 2 (Midspan) PSD comparison at 71.1 °C (160 °F) 

(Af=.25 Hz). 



Figure 5.10- Beam 2 Run 8, position 2 (midspan) PSD comparison at 93.3 °C (200°F) 

(Af=.25 Hz). 
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Figure 5.11- Beam 2 Run 8, position 2 (midspan) PSD comparison at 121.1 °C (250 °F). 



Figure 5.12- Beam 3 Run 8, position 2 (midspan) PSD comparison at 21.7 °C (71 °F) 

(Af=.25 Hz). 
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Figure 5.13- Beam 3 Run 8, position 2 (midspan) PSD comparison at 71.1 °C (160 °F) 

(Af=.25 Hz). 



Figure 5.14- Beam 3 Run 8, position 2 (midspan) PSD comparison at 93.3 °C (200 °F) 

(Af=.25 Hz). 
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Figure 5.15- Beam 3 Run 8, position 2 (midspan) PSD comparison at 121.1 °C (250 °F) 

(Af=.25 Hz). 
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Figure 5.16- Beam 2 Run 8, position 1 PSD comparison for selected temperatures 

(Af=.25 Hz). 
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Figure 5.17- Beam 3 Run 3, position 3 PSD’s comparison for selected temperatures 

(Af=.25 Hz). 



Figure 5.18- Beam 3 Run 1, position 3 PSD’s comparison through the buckling range 

(T=71°F). 
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Figure 5.19- Beam 3 Run 1, position 3 PSD’s comparison through the buckling range 

(T=100 °F). 



Figure 5.20- Beam 3 Run 1, position 3 PSD’s comparison through the buckling range 

(T=130 °F). 
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Figure 5.21- Beam 3 Run 1, position 3 PSD’s comparison through the buckling range 

(T=140 °F). 



Figure 5.22- Beam 3 Run 1, position 3 PSD’s comparison through the buckling range 

(T=150 °F). 
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Figure 5.23- Beam 3 Run 1, position 3 PSD’s comparison through the buckling range 

(T=160 °F). 
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APPENDIX A 

Projection Moire’ Interferrometry Thennoelastic Static Response Results 
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Figure 10.1- PMI full field imaging of the SMAHC Beam 2 Run 1 
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Figure 10.2- PMI full field imaging of the SMAHC Beam 2 Run 2 
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Figure 10.3- PMI full field imaging of the SMAHC Beam 2 Run 3 
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Figure 10.4- PMI full field imaging of the SMAHC Beam 2 Run 4 
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Figure 10.5- PMI full field imaging of the SMAHC Beam 2 Run 5 



5 10 

Beam Length (in) 


Figure 10.6- PMI full field imaging of the SMAHC Beam 2 Run 6 
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Figure 10.7- PMI full field imaging of the SMAHC Beam 2 Run 7. 
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Figure 10.8- PMI full field imaging of the SMAHC Beam 2 Run 8. 
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Figure 10.9- PMI full field imaging of the SMAHC Beam 2 Run 9. 
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Figure 10.10- PMI full field imaging of the SMAHC Beam 3 Run 1. 
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Figure 10.11- PMI full field imaging of the SMAHC Beam 3 Run 2 
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Figure 10.12- PMI full field imaging of the SMAHC Beam 3 Run 3 
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Figure 10.13- PMI full field imaging of the SMAHC Beam 3 Run 4 
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Figure 10.14- PMI full field imaging of the SMAHC Beam 3 Run 5 
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Figure 10.15- PMI full field imaging of the SMAHC Beam 3 Run 6 
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Figure 10.16- PMI full field imaging of the SMAHC Beam 3 Run 7 
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Figure 10.17- PMI full field imaging of the SMAHC Beam 3 Run 8. 
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APPENDIX B 

Dynamic Acceleration Data Post Processing MATLAB code 
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%This code converts an acceleration time history data file into 
%displacement time history data file by performing an FFT and then 
%integrating to get displacement and transforming back to time history. 
%Note that the code also deletes the low frequency noise present in our 
%specific data file. 

%Written by Brian Davis 

O, 

o 

clear; 

%loads data file, generates acceleration and time vectors, and defines 
%important constants 

dataname=input (' Enter the name of the data file: ', 's'); 

data= [dataname, ' . dat ' ] ; 
fid=fopen (data, ' r ' ) ; 

%f ilname=f scanf ( f id, ' %s ' , 1 ) ; 

%variables=f scanf (fid, ' %s ' , 2 ) ; 

[n, gh] =f scanf (fid, ' %g %g ' , [2 , inf ] ) ; 

%fmin=input ( ' 0 to ? frequency value is the noise present in the data?'); 

fmin=10; 

dt=n (1,2) ; 

fnyq=l / (2*dt) ; 

t=n ( 1 , : ) ' ; 

a=n (2, : ) ' ; 

N=length (a) ; 
ahat=f ft (a) ; 

%generates frequency vector 

kmin=fix (fmin*N*dt) ; 

knyq=fix (fnyq*N*dt) ; 

kl=0 : knyq; 

k2= (knyq) : — 1 : 0 ; 

k= [ kl , k2 ] ' ; 

f= (1/ (N*dt) ) . *k; 

%generates displacement coefficients by integrating, and "zeros-out" the 
%low frequency values. 

dhat= (1 . /- (2*pi . *f ) , A 2) ,*ahat; 
for j = 1 : (kmin+1) 
dhat ( j , : ) =0 ; 

end 

for j = (N- ( kmin+1 ) ) : N 
dhat ( j , : ) =0 ; 

end 

%inverse fft generates the displacement values 
x=if ft (dhat) ; 


%Compute PSD and RMS values from the generated displacement values 

fs=l/dt; 

df=l / (4096*dt) ; 

[xspec, freq] =psd(x, 4096, fs, hanning (4096) , 0) ; 
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xpsd=xspec/ fnyq; 
xrms=sqrt ( sum (xpsd) *df ) ; 


O, 

o 


%injects displacment and time vector into a matrix which is exported to a 
%new data file, note that the ' w' writes to data file and the ' t ' keeps in 
% text format so the carriage returns aren't deleted 

data_f ile= [dataname ' _disp . dat ' ] ; 
data2= [t x] ' ; 

fid=fopen (data_f ile, ' wt ' ) ; 

%fprintf ( fid, ' variables = "t (sec)" "V (volts)" \n'); 
fprintf ( f id, ' %g %g\n ' , data2 ) ; 
fclose (fid) ; 

%Write the displacement fft coefficents to data file form 

data_f ile= [dataname ' _f f t . dat ' ] ; 
data2=[f dhat] '; 
fid=fopen (data_f ile, 'wt'); 

%fprintf ( fid, ' variables = "t (sec)" "V (volts)" \n'); 
fprintf ( fid, ' %g %g\n ' , data2 ) ; 
fclose (fid) ; 

%Write the displacement psd values to data file form 

data_f ile= [dataname ' _psd . dat ' ] ; 
data2=[freq xpsd]'; 
fid=fopen (data_f ile, 'wt'); 

%fprintf ( fid, ' variables = "t (sec)" "V (volts)" \n'); 
fprintf ( fid, ' %g %g\n ' , data2 ) ; 
fclose (fid) ; 

%Write the displacement rms values to data file form 

data_f ile= [dataname ' _rms . dat ' ] ; 

data2= [xrms ] ' ; 

fid=fopen (data_f ile, 'wt'); 

%fprintf ( fid, ' variables = "t (sec)" "V (volts)" \n'); 
fprintf ( fid, ' %g %g\n ' , data2 ) ; 
fclose (fid) ; 

o, 

o 
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APPENDIX C 

Initial Imperfection Data Linearly Interpolated to Node Set, MATLAB code. 
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%Brian Davis 

%M-File that interprets experimental PMI data (deflection values in the z 
%direction) to a defined node set. 

o, 

o 


%loads experimental data file 

o, 

o 


dataname=input (' Enter the name of the dataf ile : ' , ' s ' ) ; 
data= [dataname, ' . dat ' ] ; 
fid=fopen (data, ' r ' ) ; 

[n,gh]=fscanf(fid, ' %g %g %g %g ' , [2 , inf ] ) ; 
x__data=n ( 1 , : ) ' ; 
w_data=n (2 , : ) ' ; 

o, 

o 

%generate node set 

o, 

o 

nodes= (1:185) ' ; 
for i=l : 185 

if i <= 37 

x ( i ) = ( i - 1 ) * . 5 ; 
y (i) =0; 

elseif i > 37 & i <= 74 
x ( i ) = ( i - 3 8 ) *.5; 
y ( i ) = . 2 7 5 ; 

elseif i > 74 & i <= 111 
x ( i ) = ( i - 7 5 ) *.5; 
y ( i ) = ■ 5 ; 

elseif i > 111 & i <= 148 
x (i) = (i-112 ) *.5; 
y ( i ) = . 7 2 5 ; 

else 

x (i) = (i-149) *.5; 
y (i) =1; 

end 

end 

mdl= [nodes , x ' , y ' ] ; 

O, 

o 

% interprets the data file values to a matrix corresponding to the defined 
% node sets initial deflection values (u,v,w) . 

o, 

o 

w=interpl (x_data, w_data, mdl ( : , 2 ) ) ; 
uvw= [nodes, mdl ( : , 2 : 3) . *0, w] ; 

o, 

o 

% exports to a new data file representing the initial imperfection of the 
% defined node set. 

o, 

o 

data_f ile= [dataname ' _nodeset . dat ' ] ; 
data2= [uvw] ' ; 

fid=fopen (data_f ile, ' wt ' ) ; 

%fprintf ( f id, 'variables = "node" "x coordinate" "y coordinate" "interp z 

coordinate (imperfection) \n'); 

fprintf ( f id, ' %g, %f, %f, %f \n ' , data2 ) ; 

fclose(fid); 
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APPENDIX D 

ABAQUS Static and Dynamic Analysis Input File 
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^HEADING 

Thermal post-buckling/random response of a SMAHC beam with initial 

imperfections 
* * 

*NODE, NSET=ALL_NODES 


1, 

o.. 

o.. 

0 

2, 

0.5, 

o.. 

0 

3, 

1., 

o.. 

0 

4, 

1.5, 

o.. 

0 

5, 

2., 

o.. 

0 

6, 

2.5, 

0., 

0 

7, 

3., 

0., 

0 

8, 

3.5, 

0., 

0 

9, 

4., 

0., 

0 

10, 

4.5, 

0., 

0 

11, 

5., 

0., 

0 

12, 

5.5, 

0., 

0 

13, 

6., 

0., 

0 

14, 

6.5, 

0., 

0 

15, 

7., 

0., 

0 

16, 

7.5, 

0., 

0 

17, 

8., 

0., 

0 

18, 

8.5, 

0., 

0 

19, 

9., 

0., 

0 

20, 

9.5, 

0., 

0 

21, 

10., 

0., 

0 

22, 

10.5, 

0., 

0 

23, 

11., 

0., 

0 

24, 

11.5, 

0., 

0 

25, 

12., 

0., 

0 

26, 

12.5, 

0., 

0 

27, 

13., 

0., 

0 

28, 

13.5, 

0., 

0 

29, 

14., 

0., 

0 

30, 

14.5, 

0., 

0 

31, 

15., 

0., 

0 

32, 

15.5, 

0., 

0 

33, 

16., 

0., 

0 

34, 

16.5, 

0., 

0 

35, 

17., 

0., 

0 

36, 

17.5, 

0., 

0 

37, 

18., 

0., 

0 

38, 

o.. 

0.275, 

0 

39, 

0.5, 

0.275, 

0 

40, 

1., 

0.275, 

0 

41, 

1.5, 

0.275, 

0 

42, 

2., 

0.275, 

0 

43, 

2.5, 

0.275, 

0 

44, 

3., 

0.275, 

0 

45, 

3.5, 

0.275, 

0 

46, 

4., 

0.275, 

0 

47, 

4.5, 

0.275, 

0 

48, 

5., 

0.275, 

0 

49, 

5.5, 

0.275, 

0 

o 

LO 

6., 

0.275, 

0 

51, 

6.5, 

0.275, 

0 

52, 

7., 

0.275, 

0 
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53, 

7.5, 

0.275, 

0 

■vT 

LO 

8., 

0.275, 

0 

55, 

8.5, 

0.275, 

0 

CD 

LO 

9., 

0.275, 

0 

l> 

LO 

9.5, 

0.275, 

0 

CO 

LO 

10., 

0.275, 

0 

O') 

LO 

10.5, 

0.275, 

0 

60, 

11., 

0.275, 

0 

61, 

11.5, 

0.275, 

0 

62, 

12., 

0.275, 

0 

63, 

12.5, 

0.275, 

0 

64, 

13., 

0.275, 

0 

65, 

13.5, 

0.275, 

0 

66, 

14., 

0.275, 

0 

67, 

14.5, 

0.275, 

0 

CO 

CD 

15., 

0.275, 

0 

69, 

15.5, 

0.275, 

0 

70, 

16., 

0.275, 

0 

i — 1 
[> 

16.5, 

0.275, 

0 

72, 

17., 

0.275, 

0 

73, 

17.5, 

0.275, 

0 

74, 

18., 

0.275, 

0 

75, 

o.. 

0.5, 

0 

76, 

0.5, 

0.5, 

0 

77, 

1., 

0.5, 

0 

CO 

l> 

1.5, 

0.5, 

0 

79, 

2., 

0.5, 

0 

o 

CO 

2.5, 

0.5, 

0 

1 — 1 
CO 

3., 

0.5, 

0 

CM 

CO 

3.5, 

0.5, 

0 

83, 

4., 

0.5, 

0 

CO 

4.5, 

0.5, 

0 

LO 

CO 

5., 

0.5, 

0 

CD 

CO 

5.5, 

0.5, 

0 

l> 

CO 

6., 

0.5, 

0 

CO 

CO 

6.5, 

0.5, 

0 

O'! 

CO 

7., 

0.5, 

0 

90, 

7.5, 

0.5, 

0 

1 — 1 

O') 

8., 

0.5, 

0 

92, 

8.5, 

0.5, 

0 

93, 

9., 

0.5, 

0 

94, 

9.5, 

0.5, 

0 

95, 

10., 

0.5, 

0 

96, 

LO 

O 
\ — 1 

0.5, 

0 

97, 

11., 

0.5, 

0 

CO 

O') 

11.5, 

0.5, 

0 

99, 

12., 

0.5, 

0 

100, 

12.5, 

0.5, 

0 

101, 

13., 

0.5, 

0 

102, 

13.5, 

0.5, 

0 

103, 

14., 

0.5, 

0 

104, 

14.5, 

0.5, 

0 

105, 

15., 

0.5, 

0 

106, 

15.5, 

0.5, 

0 

107, 

16., 

0.5, 

0 

108, 

16.5, 

0.5, 

0 

109, 

17., 

0.5, 

0 
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110, 

17.5, 

LO 

O 

0 

111, 

18., 

LO 

o 

0 

112, 

o.. 

0.725, 

0 

113, 

LO 

O 

0.725, 

0 

114, 

1., 

0.725, 

0 

115, 

LO 
\ — 1 

0.725, 

0 

116, 

2., 

0.725, 

0 

117, 

LO 

CM 

0.725, 

0 

118, 

3., 

0.725, 

0 

119, 

3.5, 

0.725, 

0 

120, 

4., 

0.725, 

0 

121, 

LO 

■vT 

0.725, 

0 

122, 

5., 

0.725, 

0 

123, 

LO 

LO 

0.725, 

0 

124, 

6., 

0.725, 

0 

125, 

LO 

LD 

0.725, 

0 

126, 

7., 

0.725, 

0 

127, 

LO 

[> 

0.725, 

0 

128, 

8., 

0.725, 

0 

129, 

LO 

CO 

0.725, 

0 

130, 

9., 

0.725, 

0 

131, 

9.5, 

0.725, 

0 

132, 

10., 

0.725, 

0 

133, 

LO 

o 
\ — 1 

0.725, 

0 

134, 

\ — 1 
\ — 1 

0.725, 

0 

135, 

LO 

\ — 1 
\ — 1 

0.725, 

0 

136, 

12., 

0.725, 

0 

137, 

LO 

CM 
\ — 1 

0.725, 

0 

138, 

13., 

0.725, 

0 

139, 

13.5, 

0.725, 

0 

140, 

\ — 1 

0.725, 

0 

141, 

LO 
\ — 1 

0.725, 

0 

142, 

15., 

0.725, 

0 

143, 

15.5, 

0.725, 

0 

144, 

16., 

0.725, 

0 

145, 

16.5, 

0.725, 

0 

146, 

17., 

0.725, 

0 

147, 

17.5, 

0.725, 

0 

148, 

CO 
\ — 1 

0.725, 

0 

149, 

0., 

1., 

0 

150, 

LO 

o 

1., 

0 

151, 

1., 

1., 

0 

152, 

LO 
\ — 1 

1., 

0 

153, 

2., 

1., 

0 

154, 

2.5, 

1., 

0 

155, 

3., 

1., 

0 

156, 

3.5, 

1., 

0 

157, 

4., 

1., 

0 

158, 

LO 

1., 

0 

159, 

5., 

1., 

0 

160, 

5.5, 

1., 

0 

161, 

6., 

1., 

0 

162, 

LO 

LD 

1., 

0 

163, 

7., 

1., 

0 

164, 

LO 

1., 

0 

165, 

8., 

1., 

0 

166, 

LO 

CO 

1., 

0 
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167, 

9., 

1., 

0 

168, 

LO 

CT> 

1., 

0 

169, 

10., 

1., 

0 

170, 

LO 

o 
\ — 1 

1., 

0 

171, 

11., 

1., 

0 

172, 

LO 

\ — 1 
\ — 1 

1., 

0 

173, 

12., 

1., 

0 

174, 

LO 

CM 
\ — 1 

1., 

0 

175, 

13., 

1., 

0 

176, 

13.5, 

1., 

0 

177, 

■vT 
\ — 1 

1., 

0 

178, 

LO 
\ — 1 

1-, 

0 

179, 

15., 

1., 

0 

180, 

15.5, 

1., 

0 

181, 

16., 

1., 

0 

182, 

16.5, 

1., 

0 

183, 

17., 

1., 

0 

184, 

17.5, 

1., 

0 

185, 

CO 
\ — 1 

1., 

0 


* * 

*NSET, NSET=CENTER 
93, 

*NSET, NSET=ACCELS 
84, 93, 104 
* * 

** Read geometric imperfection data from external file 
* * 

*IMPERFECTION, INPUT=beam2 run7 pmi imp nodeset.dat 
* * 


** Include temp, dependent material properties 

* * 


** GLEP 

** Date: 13-Nov-03 Time: 16:50:42 

* * 

* MATERIAL, NAME=GLEP 
*ELASTIC, TYPE=LAMINA 


** El, 

E2 , 

nul2 , 

G12 , 

G13, 

G2 3, 

T 

7 . 15e+6, 

2 . 90e+6, 

0.29, 

1 . 40e+6, 

1 . 40e+6, 

1 . 40e+6, 

60. 

7 . 15e+6, 

2 . 90e+6, 

0.29, 

1 . 40e+6, 

1 . 40e+6, 

1 . 40e+6, 

70. 

7 . 15e+6, 

2 . 90e+6, 

0.29, 

1 . 40e+6, 

1 . 40e+6, 

1 . 40e+6, 

80. 

7 . 13e+6, 

2 . 82e+6, 

0.29, 

1 . 34e+6, 

1 . 34e+6, 

1 . 34e+6, 

100 

7 . lle+6. 

2 . 75e+6, 

0.29, 

1 .29e+6, 

1 .29e+6, 

1 .29e+6, 

120 

7 . 08e+6, 

2 . 68e+6, 

0.29, 

1 ,24e+6. 

1 ,24e+6. 

1 ,24e+6. 

140 

7 . 07e+6, 

2 . 64e+6, 

0.29, 

1 . 22e+6, 

1 . 22e+6, 

1 . 22e+6, 

150 

7 . 07e+6, 

2 . 58e+6, 

0.29, 

1 . 20e+6, 

1 . 2 0e+ 6 , 

1 . 2 0e+ 6 , 

160 

7 . 06e+6, 

2 . 47e+6, 

0.29, 

1 . 15e+6, 

1 . 15e+6, 

1 . 15e+6, 

180 

7 . 05e+6, 

2 . 35e+6, 

0.29, 

1 . 10e+6, 

1 . 10e+6, 

1 . 10e+6, 

200 

7 . 05e+6, 

2 . 22e+6, 

0.29, 

0 . 98e+6, 

0 . 98e+6, 

0 . 98e+6, 

220 

7 . 04e+6, 

2 . 09e+6, 

O'! 

CM 

O 

0 . 87e+6, 

0 . 87e+6, 

0 . 87e+6, 

240 

7 . 04e+6, 

2 . 03e+6, 

0.29, 

0 . 81e+6, 

0 . 81e+6, 

0 . 81e+6, 

250 

7 . 05e+6, 

1 . 95e+6, 

0.29, 

0 . 75e+6, 

0 . 75e+6, 

0 . 75e+6, 

260 

7 . 06e+6, 

1 . 80e+6, 

0.29, 

0 . 62e+6, 

0 . 62e+6, 

0 . 62e+6, 

280 

7 . 08e+6, 

1 . 65e+6, 

0.29, 

0 . 50e+6, 

0 . 50e+6, 

0 . 50e+6, 

300 


*DENSITY 

0.00019, 
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^EXPANSION, TYPE=ORTHO, ZERO=71. 


* al. 

a2, a3, T 

0 . 2 91 9E-05 , 

0 . 6095E-05, 

60 

0 . 2 968E-05 , 

0 . 6334E-05 , 

70 

0 . 3138E-05, 

0 . 7 1 69E-05 , 

80 

0 . 3455E-05, 

0 . 9097E-05, 

100 

0 . 3664E-05, 

0 . 1060E-04 , 

120 

0 . 3750E-05, 

0 . 1149E-04, 

140 

0.37 62E-05 , 

0 . 1177E-04 , 

150 

0 . 3758E-05, 

0 . 1 1 96E-04 , 

160 

0 . 3728E-05, 

0 . 1215E-04, 

180 

0 . 3691E-05, 

0.121 9E-04 , 

200 

0 . 3666E-05, 

0 . 1220E-04 , 

220 

0 . 3665E-05, 

0 . 1222E-04 , 

240 

0 . 3673E-05, 

0 . 1225E-04 , 

250 

0 . 3687E-05, 

0 . 1230E-04 , 

260 

0 . 3723E-05, 

0 . 1244E-04 , 

280 

0 . 3763E-05, 

0 . 12 61E-04 , 

300 


** NITI 

** Date: 13-Nov-03 Time: 16:50:42 

* * 

* MATERIAL, NAME=NITI 
*ELASTIC, TYPE=LAMINA 


** El, 

E2 , 

nul2 , 

G12 , 

G13, 

G2 3, 

T 

3. 94000e+6, 

3 . 94000e+6, 

0.3, 

1 . 51535e+6, 

1 . 51535e+6, 

1 . 51535e+6, 

60. 

3. 94000e+6, 

3 . 94000e+6, 

0.3, 

1 . 51535e+6, 

1 . 51535e+6, 

1 . 51535e+6, 

70. 

3 . 59667e+6, 

3 . 59667e+6, 

0.3, 

1 . 38334e+6, 

1 . 38334e+6, 

1 . 38334e+6, 

80. 

3 . 25333e+6, 

3 . 25333e+6, 

0.3, 

1 .25128e+6, 

1 .25128e+6, 

1 .25128e+6, 

90. 

2 . 91000e+6, 

2 . 91000e+6, 

0.3, 

1 . 11923e+6, 

1 . 11923e+6, 

1 . 11923e+6, 

100 

3 . 72800e+6, 

3 . 72800e+6, 

0.3, 

1 . 43385e+6, 

1 . 43385e+6, 

1 . 43385e+6, 

110 

4 . 54600e+6, 

4 . 54600e+6, 

0.3, 

1 . 74846e+6, 

1 . 74846e+6, 

1 . 74846e+6, 

120 

5 . 36400e+6, 

5 . 36400e+6, 

0.3, 

2 . 06308e+6, 

2 . 06308e+6, 

2 . 06308e+6, 

130 

6. 18200e+6, 

6. 18200e+6, 

0.3, 

2 . 377 69e+6, 

2 . 377 69e+6, 

2 . 377 69e+6, 

140 

7 . 00000e+6, 

7 . 00000e+6, 

0.3, 

2 . 69231e+6, 

2 . 69231e+6, 

2 . 69231e+6, 

150 

7 . 95600e+6, 

7 . 95600e+6, 

0.3, 

3 . 06000e+6, 

3 . 06000e+6, 

3 . 06000e+6, 

160 

8. 91200e+6, 

8 . 91200e+6, 

0.3, 

3 . 427 69e+6, 

3 . 427 69e+6, 

3 . 427 69e+6, 

170 

9 . 31200e+6, 

9 . 31200e+6, 

0.3, 

3 . 58154e+6, 

3 . 58154e+6, 

3 . 58154e+6, 

180 

9 . 15600e+6, 

9 . 15600e+6, 

0.3, 

3 . 52154e+6, 

3 . 52154e+6, 

3 . 52154e+6, 

190 

9 . 00000e+6, 

9 . 00000e+6, 

0.3, 

3 . 46154e+6, 

3 . 46154e+6, 

3 . 46154e+6, 

200 

9 . 27000e+6, 

9 . 27000e+6, 

0.3, 

3 . 56538e+6, 

3 . 56538e+6, 

3 . 56538e+6, 

210 

9 . 54000e+6, 

9 . 54000e+6, 

0.3, 

3 . 66923e+6, 

3 . 66923e+6, 

3 . 66923e+6, 

220 

9. 81000e+6, 

9. 81000e+6, 

0.3, 

3 . 77308e+6, 

3 . 77308e+6, 

3 . 77308e+6, 

230 

10 . 0800e+6, 

10 . 0800e+6, 

0.3, 

3 . 87 692e+6, 

3 . 87 692e+6, 

3 . 87 692e+6, 

240 

10 . 3500e+6, 

10 . 3500e+6, 

0.3, 

3 . 98077e+6, 

3 . 98077e+6, 

3 . 98077e+6, 

250 

10 . 2740e+6, 

10 . 2740e+6, 

0.3, 

3 . 95154e+6, 

3 . 95154e+6, 

3 . 95154e+6, 

260 

10 . 1980e+6, 

10 . 1980e+6, 

0.3, 

3 . 92231e+6, 

3 . 92231e+6, 

3 . 92231e+6, 

270 

10 . 1220e+6, 

10 . 1220e+6, 

0.3, 

3 . 89308e+6, 

3 . 89308e+6, 

3 . 89308e+6, 

280 

1 0 . 04 60e+6 , 

1 0 . 04 60e+6 , 

0.3, 

3 . 86385e+6, 

3 . 86385e+6, 

3 . 86385e+6, 

290 

9. 97000e+6, 

9. 97000e+6, 

0.3, 

3 . 83462e+6, 

3 . 83462e+6, 

3 . 83462e+6, 

300 

9. 97000e+6, 

9. 97000e+6, 

0.3, 

3 . 83462e+6, 

3 . 83462e+6, 

3 . 83462e+6, 

320 

*DENSITY 







0.0005349, 







*EXPANSION, 

TYPE=ORTHO, 

ZERO 

=71 . 




** al 

a2 

a3 

T 




0 . 3670E-05, 

0 . 3670E-05, , 

60. 




0 . 3670E-05, 

0 . 3670E-05, , 

70. 
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-0.1430E- 

04, 

0.3670E-05, 

o 

CO 

-0.2016E- 

04, 

0 . 3670E-05, 

90. 

-0.2692E- 

04, 

0 . 3670E-05, 

100. 

-0.3082E- 

04, 

0 . 3670E-05, 

110. 

-0.5461E- 

04, 

0.3708E-05, 

120. 

-0.7336E- 

04, 

0.37 94E-05 , 

130. 

-0.7746E- 

04, 

0 . 3908E-05, 

140. 

-0.7345E- 

04, 

0 . 4 037E-05 , 

150. 

-0.6568E- 

04, 

0 . 4177E-05, 

160. 

-0.5786E- 

04, 

0 . 4324E-05, 

170. 

-0.5466E- 

04, 

0 . 4476E-05, 

180. 

-0.5452E- 

04, 

0 . 4617E-05, 

190. 

-0.5356E- 

04, 

0 . 4735E-05, 

200. 

-0.5033E- 

04, 

0 . 4836E-05, 

210. 

-0.4667E- 

04, 

0 . 4923E-05, 

220. 

-0.4348E- 

04, 

0 . 4 999E-05 , 

230. 

-0.4051E- 

04, 

0 . 5066E-05, 

240. 

-0.3765E- 

04, 

0 . 5126E-05, 

250. 

-0.3634E- 

04, 

0 . 5179E-05, 

260. 

-0.3496E- 

04, 

0.5226E-05, 

270. 

-0.3375E- 

04, 

0 . 52 69E-05 , 

280. 

-0.3257E- 

04, 

0 . 5308E-05, 

290. 

-0.3159E- 

04, 

0 . 5343E-05, 

300. 

-0.3159E- 

04, 

0 . 5343E-05, 

320. 

* * 








** SHELL 

elements (with directional properties) -- 





* * 




*SHELL SECTION, COMPOSITE, 

ELSET=SMAHCELEM 

** total 

thickness = 0.0825 

0.004875, 

3, 

GLEP, 45. 


0.006, 3, 

NITI, 0. 


0.004875, 

3, 

GLEP, -45. 


0.004875, 

3, 

GLEP, 90. 


0.004875, 

3, 

GLEP, 45. 


0.006, 3, 

NITI, 0. 


0.004875, 

3, 

GLEP, -45. 


0.004875, 

3, 

GLEP, 90. 


0.004875, 

3, 

GLEP, 90. 


0.004875, 

3, 

GLEP, -45. 


0.006, 3, 

NITI, 0. 


0.004875, 

3, 

GLEP, 45. 


0.004875, 

3, 

GLEP, 90. 


0.004875, 

3, 

GLEP, -45. 


0.006, 3, 

NITI, 0. 


0.004875, 

3, 

GLEP, 45. 



* * 


^ELEMENT, TYPE=S4, ELSET=SMAHCELEM 

* * 


37, 

38, 

39, 

76, 

75 

38, 

39, 

40, 

77, 

76 

39, 

40, 

41, 

78, 

77 

40, 

41, 

42, 

79, 

78 

41, 

42, 

43, 

80, 

79 

42, 

43, 

44, 

81, 

80 

43, 

44, 

45, 

82, 

81 
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44, 

45, 

46, 

83, 

82 

LO 

■vT 

46, 

47, 

84, 

83 

46, 

47, 

48, 

85, 

84 

47, 

48, 

49, 

86, 

85 

CO 

■vT 

49, 

50, 

87, 

86 

49, 

50, 

51, 

88, 

87 

o 

LO 

51, 

52, 

89, 

88 

\ — 1 
LO 

52, 

53, 

90, 

89 

52, 

53, 

54, 

91, 

90 

53, 

54, 

55, 

92, 

91 

54, 

55, 

56, 

93, 

92 

55, 

56, 

57, 

94, 

93 

56, 

57, 

58, 

95, 

94 

57, 

58, 

59, 

96, 

95 

CO 

LO 

59, 

60, 

97, 

96 

59, 

60, 

61, 

98, 

97 

60, 

61, 

62, 

99, 

98 

\ — 1 
LO 

62, 

63, 

100, 

99 

62, 

63, 

64, 

101, 

100 

63, 
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* * 


*SHELL SECTION, COMPOSITE, ELSET=GLEPELEM 
** total thickness = 0.078 
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3, 
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LO 
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3, 
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3, 
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3, 
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LO 

1 

0.004875, 

3, 
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0 . 

0.004875, 
* * 

3, 

GLEP, 

45. 

*ELEMENT, 

TYPE=S4 , 

ELSET=GLEPELEM 


* * 
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114, 
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152, 
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112, 

115, 

116, 

153, 
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113, 

116, 

117, 

154, 
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117, 

118, 

155, 
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115, 

118, 

119, 

156, 
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119, 
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157, 
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120, 

121, 

158, 
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118, 

121, 
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159, 
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119, 
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123, 

124, 

161, 
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121, 

124, 

125, 

162, 
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122, 

125, 

126, 

163, 
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123, 

126, 

127, 

164, 
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124, 

127, 

128, 

165, 
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125, 

128, 

129, 

166, 
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126, 

129, 

130, 

167, 

166 

127, 

130, 

131, 

168, 
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128, 

131, 

132, 

169, 
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129, 

132, 

133, 

170, 

169 

130, 

133, 

134, 

171, 
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131, 

134, 

135, 

172, 
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132, 

135, 

136, 

173, 
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133, 

136, 

137, 

174, 
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134, 

137, 

138, 
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135, 

138, 

139, 

176, 
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136, 

139, 

140, 

177, 
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137, 

140, 

141, 

178, 
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138, 

141, 

142, 
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139, 

142, 

143, 

180, 
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140, 

143, 

144, 

181, 
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141, 

144, 

145, 

182, 
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142, 

145, 

146, 

183, 
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143, 

146, 

147, 

184, 
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144, 

147, 

148, 

185, 
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* * 


*ELSET, ELSET=ALL_ELEMENTS 
SMAHCELEM, 

GLEPELEM, 

* * 

*INITIAL CONDITIONS, TYPE=TEMPERATURE , 
INPUT=beam2 run7 thermal t71 nodeset.dat 
* * 

* * 

*BOUNDARY 

CFCF, 1,6, 0. 

* * 
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* * 


*NSET, NSET=CFCF 

1, 37, 38, 74, 75, 111, 112, 148, 

149, 185 

* * 

*PSD-DEFINITION, NAME=S 1 0TO4 00 , G=1.0, INPUT=input_accel_psd . abq, 
TYPE=BASE 


** STEP data for Eigensolution at 71 F 

* * 

*STEP, PERTURBATION 

Eigensolution for SMAHC beam at 71 F 
* * 

*FREQUENCY, EIGENSOLVER=LANCZOS , NORMALIZATION=MASS 
10, , 1500. 

* * 

** Eigensolution output 
* * 

***OUTPUT, FIELD 
***NODE OUTPUT 
**U, 

*END STEP 


** STEP data for Linear Random Response at 71 F 

* * 

*STEP, PERTURBATION 

Linear random response of SMAHC beam at 71 F 
* * 

** 0-400 Hz with 400 points between eigenf requencies , 1.0 bias and linear 
frequency 
* * 

* RANDOM RESPONSE 

0. , 400., 400, 1.0, 1 
* * 

** This is only needed in the event that a subset of the modes 
** from the *FREQUENCY card are used 
* * 

*SELECT E I GENMODE S , GENERATE 

1 , 10 
* * 

*MODAL DAMPING, MODAL=DIRECT 
1,2, .0081 
3,10, .0068 
* * 

* CORRE LAT I ON , PSD=S 1 0TO4 00 , COMPLEX=NO, TYPE=CORRELATED 

2, 1. 

* * 

*DLOAD 
* * 

*BASE MOTION, DOF=3, LOAD CASE=2, TYPE=ACCELERATION 
* * 

** Random response PSD at three accel nodes 
* * 

*OUTPUT , FIELD 

*NODE OUTPUT, NSET=ACCELS 
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U, 

* * 

*END STEP 

* * 

** STEP data for post-buckling 70-80 F 

k k 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 70-80 degrees F 
* * 

* STATIC, DIRECT 

0.03,1.0 

* * 

** input thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T80 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

* * 

** STEP data for post-buckling 80-90 F 

* * 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 80-90 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7_thermal T90_nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

k k 

** STEP data for post-buckling 90-100 F 

** 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 90-100 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T100 nodeset.dat 
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* * 


** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * ** 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3, 

* * 

*END STEP 

* * 

** STEP data for post-buckling 100-110 F 

* * 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 100-110 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T110 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

k k 

** STEP data for post-buckling 110-120 F 

** 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 110-120 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T120 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

** 

** STEP data for post-buckling 120-130 F 

** 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 120-130 degrees F 
* * 
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* STATIC, DIRECT 

0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7^thermal T130 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

k k 

** STEP data for post-buckling 130-140 F 

* * 

*STEP, NLGEOM, INC=33.3333 

Post-buckling solution 130-140 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2_run7 thermal T140 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

* * 

** STEP data for post-buckling 140-150 F 

k k 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 140-150 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T150 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 
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* * 


** STEP data for post-buckling 150-160 F 

** 

*STEP, NLGEOM, INC=33.3333 
Post-buckling solution 150-160 degrees F 
* * 

* STATIC, DIRECT 
0.03,1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2_run7 thermal T160 nodeset.dat 
* * 

** Load deflection data at center node (propagates to subsequent static 
steps ) 

* * 

*OUTPUT , HISTORY 

*NODE OUTPUT, NSET=CENTER 

U3 , 

* * 

*END STEP 

* * 

** STEP data for Eigensolution at 160 F 

* * 

*STEP, PERTURBATION 

Eigensolution for SMAHC beam at 160 F 
* * 

*FREQUENCY, EIGENSOLVER=LANCZOS , NORMALIZATION=MASS 
10, , 1500. 

* * 

** Eigensolution output 
* * 

***OUTPUT, FIELD 
***NODE OUTPUT 
**U, 

*END STEP 
* * 

* * 

** STEP data for Linear Random Response at 160 F 

k k 

*STEP, PERTURBATION 

Linear random response of SMAHC beam at 160 F 
* * 

** 0-400 Hz with 400 points between eigenf requencies , 1.0 bias and linear 
frequency 
* * 

* RANDOM RESPONSE 

0. , 400., 400, 1.0, 1 
* * 

** This is only needed in the event that a subset of the modes 
** from the *FREQUENCY card are used 
* * 

*SELECT E I GENMODE S , GENERATE 

1, 10 
* * 

*MODAL DAMPING, MODAL=DIRECT 
1,2, .0017 
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3,10, .0025 
* * ** 

* CORRE LAT I ON , PSD=S 1 0TO4 00 , COMPLEX=NO, TYPE=CORRELATED 

2, 1. 

* * 

*DLOAD 
* * 

*BASE MOTION, DOF=3, LOAD CASE=2, TYPE=ACCELERATION 
* * 

** Random response PSD at three accel nodes 
* * 

*OUTPUT , FIELD 

*NODE OUTPUT, NSET=ACCELS 

U, 

* * 

*END STEP 

* * 

** STEP data for post-buckling 160-180 F 

* * 

*STEP, NLGEOM 

Post-buckling solution 160-180 degrees F 
* * 

* STATIC, DIRECT 
0.5, 1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T180 nodeset.dat 
* * 

** Load deflection data at center node (carries forward from above) 
* * 

*END STEP 

* * 

** STEP data for post-buckling 180-200 F 

* * 

*STEP, NLGEOM 

Post-buckling solution 180-200 degrees F 
* * 

* STATIC, DIRECT 

0.5, 1.0 
* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T200 nodeset.dat 
* * 

** Load deflection data at center node (carries forward from above) 
* * 

*END STEP 

** 

** STEP data for Eigensolution at 200 F 

** 

*STEP, PERTURBATION 

Eigensolution for SMAHC beam at 200 F 
* * 

*FREQUENCY, EIGENSOLVER=LANCZOS , NORMALIZATION=MASS 
10, , 1500. 

* * 
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** Eigensolution output 
* * 

***OUTPUT, FIELD 

* * *NODE OUTPUT 
**U, 

*END STEP 

* * 

** 

** STEP data for Linear Random Response at 200 F 

* * 

*STEP, PERTURBATION 

Linear random response of SMAHC beam at 200 F 
* * 

** 0-400 Hz with 400 points between eigenf requencies , 1.0 bias and linear 
frequency 
* * 

* RANDOM RESPONSE 

0. , 400., 400, 1.0, 1 
* * 

** This is only needed in the event that a subset of the modes 
** from the *FREQUENCY card are used 
* * 

*SELECT E I GENMODE S , GENERATE 

1, 10 
* * 

*MODAL DAMPING, MODAL=DIRECT 
1,2, .0016 
3,10, .0016 
* * 

* CORRE LAT I ON , PSD=S 1 0TO4 00 , COMPLEX=NO, TYPE=CORRELATED 

2, 1. 

* * 

*DLOAD 
* * 

*BASE MOTION, DOF=3, LOAD CASE=2, TYPE=ACCELERATION 
* * 

** Random response PSD at 3 accel nodes 
* * 

*OUTPUT , FIELD 

*NODE OUTPUT, NSET=ACCELS 

U, 

* * 

*END STEP 

* * 

** STEP data for post-buckling 200-220 F 

* * 

*STEP, NLGEOM 

Post-buckling solution 200-220 degrees F 
* * 

* STATIC, DIRECT 

0.5, 1.0 

* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T220 nodeset.dat 
* * 

** Load deflection data at center node (carries forward from above) 
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* * 


*END STEP 

** 

** STEP data for post-buckling 220-250 F 

* * 

*STEP, NLGEOM 

Post-buckling solution 220-250 degrees F 
* * 

* STATIC, DIRECT 

0.5, 1.0 

* * 

** Uniform thermal load 
* * 

*TEMPERATURE, INPUT=beam2 run7 thermal T250 nodeset.dat 
* * 

** Load deflection data at center node (carries forward from above) 

* * 

*END STEP 

* * 

** STEP data for Eigensolution at 250 F 

* * 

*STEP, PERTURBATION 

Eigensolution for SMAHC beam at 250 F 
* * 

*FREQUENCY, EIGENSOLVER=LANCZOS , NORMALIZATION=MASS 

10, , 1500. 

* * 

** Eigensolution output 
* * 

* * *OUTPUT , FIELD 
***NODE OUTPUT 
**U3, 

*END STEP 

* * 

* * 

** STEP data for Linear Random Response at 250 F 

* * 

*STEP, PERTURBATION 

Linear random response of SMAHC beam at 250 F 
* * 

** 0-400 Hz with 400 points between eigenf requencies , 1.0 bias and linear 
frequency 
* * 

* RANDOM RESPONSE 

0. , 400., 400, 1.0, 1 
* * 

** This is only needed in the event that a subset of the modes 
** from the *FREQUENCY card are used 
* * 

*SELECT E I GENMODE S , GENERATE 

1, 10 
* * 

* MO DAL DAMPING, MODAL=DIRECT 
1,2, .0021 
3,10, .0022 
* * 

* CORRE LAT I ON , PSD=S 1 0TO4 00 , COMPLEX=NO, TYPE=CORRELATED 
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2 , 1 . 

* * 

*DLOAD 
* * 

*BASE MOTION, DOF=3, LOAD CASE=2, TYPE=ACCELERATION 
* * 

** Random response PSD at 3 accel nodes 
* * 

*OUTPUT , FIELD 

*NODE OUTPUT, NSET=ACCELS 

U, 

* * 

*END STEP 
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